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ABSTRACT 

This  thesis  presents  the  results  of  numerical  simulations  of  mode-locked  figure 
eight  lasers  and  their  components:  fiber  amplifiers  and  nonlinear  optical  loop  mirrors 
(NOLMs).  The  computations  were  designed  to  study  pulse  evolution  in  optical 
amplifiers  and  NOLMs  with  periodic  repetition  of  these  elements.  Since  fiber  laser 
systems  also  include  birefringent  fiber,  the  effects  of  fiber  birefringence  was 
incorporated  into  the  simulations. 

My  studies  of  pulse  amplification  in  non-birefringent  amplifiers  show  pulse 
breakup  when  their  energies  exceed  4.5  fundamental  soliton  energies.  In  birefringent 
fibers  pulse  breakup  is  also  found,  but  the  two  orthogonally  polarized  pulses  propagate 
together.  I  find  that  their  behavior  is  related  to  the  properties  of  a  vector  soliton. 

I  found  that  vector  waves  have  close  to  unity  transmission  through  a  birefringent 
NOLM,  but  the  pulse  shape  is  distorted.  This  shape  distortion  reduces  subsequent 
transmissions  through  the  NOLM.  The  energy  required  for  peak  transmission  of  the 
pulse  is  predicted  by  the  theory  based  on  vector  solitons.  The  same  theory  also 
predicted  the  low  intensity  transmission.  The  performance  of  the  NOLM  with 
birefringent  fiber  could  not  be  improved  by  altering  the  polarization  state  of  the  pulse 
from  linear  polarization;  the  polarization  controller  introduced  pulse  distortion  that 
resulted  in  excessive  loss. 

I  found  an  instability  in  the  steady- state  operation  of  the  figure  eight  laser,  which 
is  due  to  pulse  reshaping  during  propagation  in  the  amplifier  section.  To  remove  this 
instability  I  introduced  the  concept  of  dispersion  balancing;  by  increasing  the  dispersion 
in  the  amplifier  section,  the  pulse  can  propagate  nearly  as  a  fundamental  soliton  in  both 
the  amplifier  and  the  NOLM  sections  of  the  laser.  This  eliminated  a  major  source  of 
dispersive  wave  shedding  and  allowed  the  laser  operation  to  become  independent  of  the 
amplifier  length.  Sidebands  were  found  on  the  pulse  spectrum  and  their  maxima 
corresponded  well  with  the  periodic  resonance  model. 
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Chapter  One  :  Introduction 


This  thesis  describes  the  results  of  simulations  of  the  figure  eight  laser  and  its 
components.  The  intent  has  been  to  improve  the  understanding  of  the  behavior  of  these 
lasers  and  to  propose  improvements  to  their  design.  A  great  deal  of  the  behavior  of 
these  lasers  appears  to  be  controlled  by  the  polarization  state  of  the  light.  For  this 
reason,  the  evolution  of  the  polarization  state  has  been  extensively  studied.  The  laser 
simulated  is  shown  in  Figure  1.1.  It  consists  of  an  output  coupler,  an  optical  isolator,  an 
amplifier,  and  a  Nonlinear  Optical  Loop  Mirror  (NOLM)  which  serves  to  mode-lock  the 
laser.  The  thesis  examines  how  each  of  the  components  can  be  simulated,  and  then  how 
simulations  of  the  parts  are  used  to  model  the  total  system.  Propagation  of  pulses  in 
fiber  and  their  simulation  is  discussed  in  Chapter  2.  The  effects  of  birefringence  are 
introduced  in  Chapter  3.  Chapter  4  describes  how  the  amplifier  is  simulated  and  the 
NOLM  is  discussed  in  Chapter  5.  In  Chapter  6,  all  these  components  are  drawn 
together  and  simulations  of  the  laser  system  are  described. 

It  is  hoped  that  these  mode-locked  fiber  lasers  will  provide  a  convenient  way  of 
producing  short  optical  pulses  at  high  repetition  rates.  Optical  fiber  offers  an  attractive 
means  of  producing  many  useful  optical  effects.  These  properties  can  be  used  to 
produce  fiber  lasers  and  fiber  devices.  Since  the  invention  of  the  laser,  optics  has  been 
viewed  as  an  attractive  means  for  very  high  speed  data  transmission.1  This  is  due  to  the 
extremely  large  bandwidths  available,  compared  to  traditional  electronics.  Because  of 
the  poor  properties  of  optical  transmission  through  the  atmosphere  and  the  inevitable 
tendency  of  optical  beams  to  diverge  due  to  diffraction,  the  development  of  optical 
communications  had  to  wait  for  the  invention  of  low  loss  optical  fibers.  In  addition  to 
producing  a  low  loss  transmission  path,  fiber  also  produces  amplified  nonlinear  effects. 
The  fiber  continually  focuses  down  the  light,  overcoming  the  tendency  of  light  to  spread 
out  due  to  diffraction.  Apart  from  enabling  the  fiber  to  act  as  a  means  of  guiding  light 
from  one  point  to  another,  this  means  that  light  in  a  fiber  is  of  moderately  high  intensity 
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for  long  distances.  Since  the  distances  can  be  made  long,  large  nonlinear  optical  effects 
can  be  produced  in  standard  silica  fibers  in  spite  of  relatively  low  nonlinear  coefficients. 
By 

polarization 


Figure  1.1  NOLM  Figure  Eight  Laser 


comparison,  an  experiment  involving  light  in  free  space  and  bulk  nonlinear  optical 
material,  requires  the  light  to  be  focused  down  by  a  lens.3  The  light  near  the  focal  point 
of  the  lens  can  attain  very  high  intensities,  but  only  for  a  very  short  distance.  Indeed,  the 
more  the  light  is  focused  down  to  get  high  intensities,  the  shorter  this  distance 
becomes.4  These  facts  make  optical  fibers  interesting  as  a  system  for  investigating 
practical  nonlinear  optics,  involving  modest  amounts  of  optical  power.  The  same 
optical  confinement  is  needed  to  make  efficient  optical  amplifiers.  This  means  that 
fibers  can  make  excellent  optical  cavities  for  lasers. 

The  propagation  of  light  in  a  nonlinear  medium  makes  possible  a  variety  of 
useful  phenomena.  Since  the  light  can  interact  with  light  from  other  sources,  all-optical 
switches  are  possible.5  Recently,  fiber  amplifiers  have  emerged  as  a  practical  device.6,7 

The  nonlinear  effect  in  fibers  results  in  the  formation  of  optical  pulses  known  as 
optical  solitons.8  In  a  medium  in  which  there  is  anomalous  dispersion,  where  the  group 
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velocity  increases  with  frequency,  the  medium  will  support  solitons.  In  a  dispersive 
fiber,  if  one  ignores  nonlinear  effects,  a  short  optical  pulse  will  tend  to  spread  in  time 
since  light  of  different  frequencies  will  propagate  at  different  velocities.  Optical 
solitons  exist  when  the  nonlinearity  of  the  fiber  cancels  out  the  effect  of  dispersion.  In 
silica  fiber,  there  is  anomalous  dispersion  of  light  with  a  wavelength  of  1.55  pm.  This 
also  corresponds  with  the  wavelength  which  experiences  minimum  loss  in  silica  fiber. 
Since  these  pulses  remain  stable,  they  are  of  interest  for  long  distance  very  high  speed 
communications. 

Optical  nonlinearity  in  fiber  has  been  used  to  produce  a  variety  of  all  optical 
switches.5  Loop  mirrors  work  by  using  the  phase  difference  built  up  as  light  of  unequal 
intensities  propagates  in  a  Sagnac  loop.9  Dual  core  fibers  have  been  made  in  which 
light  in  the  two  cores  can  interact.  This  interaction  may  be  used  to  cause  switching 
behavior.10  Interaction  of  light  of  different  polarizations  in  birefringent  fiber  causes  the 
frequency  of  the  light  in  both  polarizations  to  shift  so  that  pulses  of  light  in  both 
polarizations  may  propagate  as  one  coupled  soliton.  This  phenomenon,  known  as 
soliton  dragging,  has  been  used  to  produce  a  switch.11  In  a  birefringent  fiber,  the 
polarization  of  elliptically  polarized  light  evolves  at  a  rate  proportional  to  its  intensity, 
so  an  optical  polarizer  can  be  placed  in  the  system  to  produce  intensity  dependent 
switching.  Time  division  multiplexing  has  been  demonstrated  using  these  devices. 

The  combination  of  optical  amplifiers  and  intensity  dependent  switching  makes 
possible  many  forms  of  mode-locked  fiber  lasers  for  the  production  of  short  optical 
pulses6.  Indeed,  lasers  using  Sagnac  switches  and  polarization  evolution  have  been 
demonstrated  experimentally. 12,13,14  Dual  core  fiber  lasers  have  been  proposed  and 
studied  by  simulations.15  The  appeal  of  an  all  fiber  laser  is  that  it  would  appear  to  offer 
a  robust  and  inexpensive  means  of  producing  short  optical  pulses.  All  these  are  also 
made  possible  by  the  large  bandwidth  of  Erbium  doped  fiber  amplifiers.  This  wide 
bandwidth  also  means  that  the  lasers  have  a  very  large  saturation  energy  and  this  means 
pulses  produced  in  fiber  lasers  tend  to  be  highly  symmetric  and  unchirpped.  These 
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advantages,  especially  as  compared  to  semiconductor  lasers,  have  led  commercial  firms 
to  examine  fiber  lasers  as  a  means  of  producing  pulses  for  long  distance 
communications.  In  this  application,  the  pulses  formed  must  be  close  to  solitons;  any 
imperfections  will  result  in  unwanted  background  noise.  Alternatively,  the  wide 
bandwidth,  in  Erbium  doped  fiber  amplifiers,  has  been  used  to  produce  tunability  over  a 
bandwidth  of  about  30  nm. 

This  thesis  will  examine  a  number  of  these  devices  by  means  of  numerical 
simulation.  In  the  second  chapter,  the  basic  physics  of  pulse  propagation  in  fiber  will  be 
examined  along  with  the  numerical  routines  used  to  simulate  it.  The  third  chapter  will 
extend  these  topics  to  the  case  of  birefringent  fibers.  In  fact  most  fibers  are  birefringent, 
so  this  topic  must  be  taken  up.  Fiber  amplifiers  are  described  and  the  model  used  in  this 
thesis  is  discussed  in  chapter  four.  The  simulations  discussed  here  extend  the 
understanding  of  pulse  amplification  in  birefringent  amplifiers.  These  chapters  lay  the 
groundwork  for  the  simulation  of  other  fiber  devices. 

Two  systems  have  been  examined  in  detail.  They  are  the  Nonlinear  Optical 
Loop  Mirror  (NOLM),  and  the  figure  eight  laser.  The  NOLM  is  a  switch  based  on  the 
Sagnac  effect  mentioned  above.  In  addition  to  its  other  uses,  it  can  be  used  as  a  means 
of  mode-locking  a  fiber  laser.  The  understanding  of  the  NOLM  has  been  extended  by 
introducing  the  effect  of  birefringence.  I  have  examined  configurations  which  produce 
maximum  transmission  along  with  the  effect  of  the  switch  on  the  pulses  which  are  being 
transmitted. 

The  simulations  of  the  fiber  laser  carried  out  for  this  thesis  extend  previous  work 
in  a  number  of  ways.  I  have  introduced  the  effect  of  propagation  in  the  fiber  amplifier. 
This  effect  has  proven  to  be  of  great  importance  to  the  understanding  of  the  laser.  These 
lasers  are  marked  by  a  number  of  different  stability  regimes.  These  appear  in  the 
simulations.  They  have  been  examined,  and  in  one  case  a  means  for  eliminating  one 
form  of  the  instability  has  been  found.  In  the  final  chapter,  the  conclusions  of  the 
research  are  drawn. 


4 


Chapter  Two  :  Simulation  of  pulse  propagation 


This  chapter  is  intended  to  serve  as  an  introduction  to  the  physics  of  pulse 
propagation  in  optical  fiber.  The  intention  is  to  show  how  the  basic  equations  are  derived 
and  to  discuss  the  behavior  of  optical  solitons.  The  major  concern  is  the  behavior  of  light 
in  single-mode  step-index  silica  fiber1.  This  is  the  sort  of  fiber  used  in  building  the 
optical  devices  described  in  the  rest  of  the  thesis.  Single  mode  fiber  offers  the  advantage 
that  transverse  modes  propagate  with  different  velocities  in  the  fiber  so  that  a  multimode 
signal  will  break  up  as  it  propagates.  A  step  index  fiber  has  a  cylindrical  core  with  an 
index  of  refraction  nu  in  which  light  is  largely  confined.  This  is  surrounded  by  a  layer  of 
silica  referred  to  as  the  cladding  which  has  an  index  of  refraction  n2.  If  n2<m,  light  will 
tend  to  be  confined  to  the  core  and  the  fiber  will  act  as  an  optical  wavguide. 


Linear  Pulse  Propagation  in  a  Fiber 

The  fundamental  component  of  a  simulation  of  fiber  lasers  is  the  simulation  of 
pulse  propagation  in  fiber.  Simulation  of  optical  pulses  is  made  more  manageable  if  the 
function  describing  the  electric  field  is  broken  down  into  components.  The  electric  field 
of  the  pulse  can  be  described  in  the  following  fashion1 

E(r,t)  =  xA(t)F(x,y)exp{i[co0t  -  Pz]} .  (2.1) 

In  this  equation,  F(x,y)  describes  the  transverse  shape  of  the  pulse  in  the  fiber.  The 
exponential  term  describes  the  rapid  changes  of  the  electric  field.  The  parameter  (00  is  the 
central  frequency  and  p  is  the  propagation  constant.  The  function  A(t)  represents  the 
envelope  function,  that  is,  the  pulse  shape.  The  functions  F(x,y)  and  p  are  set  by 
properties  of  the  fiber  such  as  the  radius  of  the  core,  the  difference  in  the  index  of 
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refraction  between  the  core  and  the  cladding,  and  by  the  central  frequency.  The  envelope 
function  is  the  part  which  evolves  during  propagation  through  fiber,  so  it  is  the  focus  of 
further  investigation. 

The  evolution  of  the  envelope  function  is  studied  using  the  nonlinear  Schrodinger 
Equation  (NLSE).  The  linear  equations  describing  propagation  in  a  fiber  are  derived  first. 
The  nonlinear  equations  are  derived  as  a  perturbation  to  the  linear  case.  This  derivation 
can  be  found  in  standard  references.  ’  The  starting  point  for  its  development  is  the  wave 
equation 


V2E(f,t) 


1  32E(f,t) 
c2  di1 


52PlM 

at2  ’ 


(2.2) 


where  PL  is  the  linear  polarization  and  P,  =  e0%E.  We  will  assume  the  electric  field  has 
the  form  in  Equation  2.1,  but  it  will  be  considered  in  frequency  space  where  it  becomes 

E(r,co)  =  F(x,y)A(co  -co0)exp[i(Pz)].  (2.3) 


If  this  is  inserted  in  Equation  2.2,  we  get 


d2F(x,y)  a2F(x,y) 
Sx2  dy2 


+  p0£ocooF(x,y)- p2F(x,y)  =  -p0s0xF(x,y)  . 


(2.4) 


Collecting  terms  this  becomes 


koF(x,y)-p2F(x,y)  =  0  . 


(2.5) 


This  equation  can  be  transformed  into  cylindrical  equations.  The  electric  field  then 
becomes 

E(p,z,©)  =  A((0)F(p)e(im*)e-ipz, 
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and  Equation  2.5  becomes 


d2F  1  dF 

— T  - 

dp  p  dp 


p2  2 


F  =  0, 


(2.6) 


where 


K2=n2k2-P2. 

The  variable  n,  is  the  refractive  index  which  takes  on  values  ni  in  the  core  and  n2  outside 
the  core.  Since  Equation  2.6  is  a  Bessel  equation,  its  solutions  will  be  various  forms  of 
the  Bessel  function.  The  solutions  are  further  constrained  by  the  facts  that  the  field  must 
be  finite  at  the  origin  and  go  to  zero  as  p  goes  to  infinity.  The  field  must  also  be 
continuous  at  the  boundary  of  the  core  and  the  cladding.  The  first  two  boundary 
conditions  require  that  within  the  core  the  solution  will  be 

p'-'fp) = j„Fp) 

inside  the  core  and 

F|n)(p)=  K„(yp) 

outside  the  core,  where  Jm  is  the  Bessel  function,  Km  is  the  modified  Bessel  function  and 

y=(^-nX)y\ 

The  second  boundary  condition  serves  to  define  an  eigenvalue  equation  that  sets  the 
values  of  (3,  y  and  k  for  a  given  set  of  values  of  the  core  radius,  ni,  n2,  A,o,  and  the  order  of 
the  solution,  m.  Values  of  the  propagation  constant  of  order  m  will  be  expressed  as  P(m). 
The  normalized  frequency,  V,  is  defined  by 
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V  =  k0a(n2  -n2)^. 

where  a  is  the  core  radius.  The  fiber  will  only  support  one  mode  if  V  is  less  than  about 
2.405. 


Nonlinear  Pulse  Propagation 


To  study  the  nonlinear  phenomena  in  a  fiber,  one  starts  with  the  nonlinear  wave 
equation 


V2E 


1  82E  52Pl  d2?^ 

c2  a2  ~  a2  ~  a2  ■ 


(2.7) 


In  the  last  equation,  PL  is  the  linear  part  of  the  polarization  and  PNI  is  the  nonlinear  part. 
We  will  transform  this  equation  by  looking  for  solutions  of  the  form  of  Equation  2.3, 

E(r,co)  =  xF(x,y)A(©  -co0,z)exp[i(P0z)] .  (2.8) 


The  variable  p0  is  the  propagation  constant.  The  envelope  function  is  now  a  function  of 
z,  as  well  as  co  because  the  nonlinear  part  of  the  polarization  can  cause  it  to  evolve  as  it 
propagates  in  the  fiber.  When  this  is  put  into  2.7  we  get 


a2F 

dx2 


A  +  -^A  +  ^41F  +  2ip0 


dy 2 


dz 


—  F  +  s((»)k2FA  -  p02FA  =  0 . 
dz 


(2.9) 


This  equation  is  separable.  If  we  use  P 2  as  the  separation  constant,  which  we  will  obtain 
below  as  the  result  of  solving  an  eigenvalue  problem,  the  result  of  the  separation  is 
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(2.10) 


32F  52F 


dx  dy 


+  ^T  +  s(co)k2F-p2F  =  0 


8  A  . D  5 A  Dj~  0  2 'j'  a 
'VT  +  21Po—  +  P  A-P0  A-0. 
5z  oz 


(2.11) 


We  will  first  examine  Equation  2.10.  This  will  give  us  the  form  of  the  transverse 
functions.  To  go  further  we  need  to  get  Equation  2.10  into  a  form  comparable  to  that  of 
2.5.  To  do  this  we  will  use  the  following  approximation 


s(co)  =  (n  +  An)2  =  n2  +  2nAn 


(2.12) 


if  An  is  small. 


The  quantity  An  combines  the  nonlinear  and  imaginary  parts  of  the  index  of 


refraction 


An  =  nkI  +  ; 


(2.13) 


where  I  is  the  optical  intensity.  It  is  assumed  that  the  two  additional  effects  that 
contribute  to  An  are  attenuation  and  the  optical  Kerr  effect.  The  first  term  of  the  last 
equation  describes  the  optical  Kerr  effect.17  The  second  describes  the  effect  of 
attenuation  or  amplification  in  the  fiber.  The  intensity  can  be  defined  as 


|F|2  A2 


r00  r°°  i  1 9 

F  dxdy 

J-00  J-00  1  1 


If  we  put  Equation  2.13  into  2.10,  we  get 


a2F  5  F 


dx  dy 


+  ^-y  +  n2k2F  +  2nAnk2F-p2F  =  0. 


(2.14) 
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Equation  2.14  is  a  perturbed  form  of  2.5.  It  can  be  solved  by  first  order  perturbation 
theory.18,19  We  seek  a  value  of  the  eigenvalue  of  the  form 


P  =  P(©)  +  XTAp,  (2.15) 

and  eigenfunctions  of  the  form 

F  =  F0+TF1.  (2.16) 

Where  F0  are  the  eigenfunctions  of  the  unperturbed  problem,  Fi  is  the  function  which  is 
the  first  correction,  and  F  is  the  perturbation  parameter.  The  perturbation  parameter  is 
proportional  to  the  nonlinear  coefficient  nk  and  the  absorption  coefficient  a,  which  are 
small  parameters.  If  these  are  substituted  into  2.14  and  the  result  is  separated  into  like 
powers  of  F,  we  get 

|SL  +  0  +  nJkSFo-P(o))!Fo=O  (2.17) 

ox  oy 

and 

52F  52F 

+  -^L  +  n2k2F1+2nAnk2F0-2p(©)ApF0-p(©)2F1  =0.  (2.18) 

Equation  2.17  is  the  linear  propagation  equation  derived  above.  The  solutions  to  this 
equation  will  be  the  same  as  before.  The  values  of  Fj  are  sought  in  terms  of  a  series 
expansion  in  terms  of  the  complete  solutions  of  2.17,  that  is  in  terms  of 

F0(m)  m=  1,2,3,...  . 

The  expansion  can  be  expressed  as 
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f, = 2><->Fr 


This  is  substituted  into  2.18  to  obtain 


-i  (m)  ^  (m) 

aF»  ^F» 


-+- 


dxl  dy1 


'-o-Lo 


+  2nAnk’F'”-[p(")(a))]JXa(")Fr,=2[pl")(ot,)]APF'-1. 


The  last  equation  reduces  to 

^a(m)[p(m)(co)]F0(m)  +2nAnk20F0(n)  -[p(n)(co)]2^a(m)F0(m>  =  2[p^(co)j  ApF0ln> . 

in  m 

This  equation  is  premultiplied  by  F0(k)  and  integrated  over  all  space.  Since  the  functions 
Fo(m>  are  orthogonal,  we  are  left  with 


,oo 


[p(k)(co)]+2nk2  f  f  [F(k)]*AnF(n)dxdy-a(k)[p(n)(co)]  =2[p(n)(co)]  ApH[F(k)]'F(n)dxdy. 


In  perturbation  analysis,  one  would  expand  F  \  as  a  power  series  of  the  eigenfunctions  of 
Fo.  This  case  is  simple  as  there  is  only  one  mode.  If  we  insert  Fo(1)  into  the  last  Equation, 
and  define  p(l)(©)=P((o),  the  result  is 


nk2  ]][F<k)]*AnF<1)dxdy  =  P(e>)  Ap  ]][F'k)]*F<1)dxdy . 


(2.19) 


Solving  for  Ap  and  dropping  the  superscripts  this  becomes 


k0  J  |F0*AnF0dxdy  k0nk|A|  J  J|F0|4dxdy 


AP: 


J  jFo’F0dxdy  j  jF0*F0dxdy 


ia 

- O - 1 - 5 

A2  2 


(2.20) 


where  we  have  made  the  assumption  that  P(co)  =  nk 


o- 
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We  may  now  use  the  eigenvalue  P  in  Equation  2.11.  We  will  also  make  the 
slowly  varying  envelope  approximation  and  eliminate  the  term  involving  the  second 
derivative.  We  will  further  use  the  approximation 


P2-P02  =(P-P0)‘-2PS+2PP0S2P0(P-P0)- 


(2.21) 


This  is  justified  since  the  difference  is  small.  This  leaves  us  with 


2ip0— +  2Po[P(co)  +  Ap-po]A  =  0, 
oz 


(2.22) 


which  can  be  expressed  as 


f5A 

— -  =  i[P(oi)  +  Ap-p0]A. 

oz 


(2.23) 


The  function  P(a>)  can  be  approximated  by  its  Taylor  expansion 


1  2  1 

p(co)  =  Po  +  P,(<0  -  (00)  +  -  p2(©  -  co0)2  +  -  P3(co  -  co0)3+. . 


(2.24) 


where 


This  is  substituted  into  2.23  to  get 


8A  1  ?  1  _ 

™  =  i[Pi(co-ro0)  +  ~P2((D-co0)  +-p3((D-co0)  +Ap]A. 


(2.25) 


The  inverse  Fourier  transform  of  this  is  taken  noting  that  the  inverse  Fourier  transform  of 


(go-cdo)  is  i  The  result  is 
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(2.26) 


dA  5A  i  d2 A  p3  a3A  ,ArJA 

dz  Pl  a  2 P2  a2  6  a3  p 


When  2.20  is  substituted  into  2.26,  we  get 


dA  „  dA  i  „  52A  p3  53A  a  a  .  ,  A  |2  A 

—  =  — B,  —  —  —  B-,  —  ,  +  —  A  +  iy  A  A, 

dz  Pl  a  2 P2  a2  6  a3  2  n 


(2.27) 


where  the  nonlinear  coefficient  y  is  given  by 


nk©o 


(2.28) 


\-oo  -co 
^  eff  —  oo.  co. 


J  J|F|2  dxdy 

-co  -co  ' 

HlFtdxdy 


(2.29) 


Some  additional  terms  can  be  added  which  include  effects  found  in  short  pulses. 

2 

In  the  first  place  a  more  extended  expansion  of  [3(g))  than  was  given  in  2.24  can  be  used. 

2 

One  may  expand  [3(o))  in  terms  of  co  and  in  terms  of  |A|  giving 


P(o>)  =  P0  +P!(co -co0)+  ip2(<0 -co0)z  +tp,(co-(00)J  + 
dA  dA  5® 


(2.30) 


If  this  is  put  into  2.23  and  the  inverse  Fourier  transform  is  taken,  we  obtain 


dA  n  dA  i  d2 A  p3  S3A  a  .  Ia,2a  3  /,  ,2  \ 

&  -  sT  -  2  +  '  2  ■ A  +  ,riA|  A  ‘  a'  a  (|Ai  A)>  (231) 
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1  J 

where  a!  =  This  new  term  is  responsible  for  the  phenomenon  know  as  self- 

steepening  which  will  be  described  later. 


An  additional  term  is  caused  by  the  Raman  effect.  This  is  phenomenon  known  as 
the  self-frequency  phase  shift.  The  Raman  effect  has  been  found  to  have  a  form5 

1 ~dz  =  A  jfrMA|2dco-  (2.32) 

—oo 


The  response  function  fr  can  be  approximated  by  the  linear  term  of  its  Taylor  expansion, 
If  the  inverse  Fourier  transform  of  the  last  expression  is  taken  we  get 


(2.33) 


The  constant  a2  is  an  experimentally  measured  value.  It  is  more  conventional  to  refer  to 
the  Raman  lifetime  a2  as  Tr. 

The  derivation  above  results  in  the  extended  version  of  the  nonlinear  Schrodinger 
equation  (NLSE)21,22  which  can  be  used  to  describe  pulse  evolution  for  pulses  as  short  as 
about  10  fs.  In  some  places,  constants  will  be  changed  to  make  the  equations  more 
useful.  It  is  given  by 


where  T  is  time,  as  measured  in  a  reference  frame  moving  along  with  the  pulse.  That  is 
z 

T  =  t  ~  —  where  t  is  time  and  vg  is  group  velocity.  The  parameters  p2  and  p  are  the 

vg  3 

second  and  third  order  dispersions  of  the  fiber. 
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It  is  often  easier  to  work  with  the  nonlinear  Schrodinger  equation  in  normalized 
units.  This  equation  is  often  given  in  normalized  units 


.  du  1  S  u  8  u 

1  —  H - r-  -  id  — -  = 

d ^  2  dx2  dx 3 

-m2  ||2  .  d  (,  a  \  d 

-N  u  u  +  is —  u  u  -tru  — 

/  fix 


(l“l2)  ■ 


(2.35) 


The  variables  U,  x,  and  are  given  by 


VpT 


where  PQ  is  the  peak  power  of  the  pulse.  The  dispersion  length  is  represented  by  Ld, 

T2 

where  LD  =  — ^ .  The  pulse  width  measured  as  the  full  width  at  half  maximum  is  given 


by  Tfwhm-  The  pulse  width  normally  used  is  T  ,  given  by  T 


1.763  T  for  a 


hyperbolic  secant  shaped  pulse.  The  pulse  will  be  significantly  distorted  by  linear 
dispersion  if  it  propagates  over  a  distance  long  relative  to  a  dispersion  length.  The 
parameters  8,  tr,  and  s  describe  the  effects  of  higher-order  dispersion,  Raman  scattering 
and  self-steepening,  respectively.  They  are  given  by 

s .  P3  =J_  2 

6|P2|T0  ’  Tr  Tr’  s  oj0T0 

The  parameter  N  is  given  by 


N2  =lZoZk. 


(2.36) 


It  should  be  noted  that  soliton  distances  are  sometimes  given  in  terms  of  the 
soliton  period,  which  is 
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(2.37) 


D  ' 


Normalized  distances  are  also  given  in  terms  of  dispersion  lengths.  This  second  system 
will  be  used  in  this  thesis,  except  when  citing  the  work  of  others  or  when  the  components 
being  simulated  had  been  described  by  others.  In  this  way  the  simulations  are  directly 
comparable. 

When  higher  order  effects  are  ignored,  Equation  2.35  reduces  to 


.  du  1  d2  u 

i - 1 - t- 

d %  2  dx2 


This  equation  is  known  to  have  soliton  solutions.1,16  The  second  term  on  the  left  hand 
side  of  the  last  equation  describes  the  effect  of  second  order  linear  dispersion.  The  term 
on  the  right  hand  side  describes  the  effect  of  the  optical  Kerr  nonlinearity  on  the  pulse 
envelope.  Optical  solitons  exist  when  these  two  effects  are  balanced.  Solitons  have  the 
property  that  the  envelope  function  evolves  in  a  periodic  fashion  as  it  propagates.4  When 
the  pulse  starts  with  a  form 

u(0,t)  =  Nsech(x) 


and  N  is  an  integer,  the  pulse  will  act  as  a  soliton  with  a  period  zo.  The  work  in  this 
thesis  focuses  on  the  fundamental  soliton,  that  is  when  N=1 .  In  practice  it  has  been  found 
-that  if  0.5  <  N  <  1 .5,  the  pulse  will  evolve  into  a  fundamental  soliton.  The  envelope 
function  of  a  fundamental  soliton  is 


u(^,t)  =  2q  sec  hurl'll)  exp(2ir|2£,). 
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It  should  be  noted  that  the  amplitude  of  the  soliton,  its  width,  and  how  its  phase  changes 
as  it  propagates  are  all  linked. 

The  work  reported  in  this  thesis  is  expressed  in  normalized  units  of  the  sort 
accompanying  Equation  2.35.  At  times  it  is  necessary  to  choose  a  standard  soliton  width, 
in  order  to  convert  to  physical  units.  The  results  have  been  normalized  to  a  300  fs 
soliton,  so  the  times  reported  here  should  be  multiplied  by  300  fs  to  be  converted  to  real 
time.  The  distances  are  usually  given  in  terms  of  the  dispersion  length  which  is  about  4.5 
m  for  the  data  expressed  above.  The  energy  of  the  soliton  will  be  26  pJ.  Other  units  can 
be  derived  from  these. 
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Table  2.1  Typical  data  for  optical  fibers 


Tr 

6  fsec 

P2 

P3 

n  core  index 

1.46 

An,  core-cladding  index  difference 

0.0036 

nk 

a,  core  radius,  standard  fiber 

4.15pm 

Additional  Data  for  Typical  Doped  Fiber 

a,  core  radius 

3.25pm 

Na 

■ 

An,  core  cladding  index  difference 

0.0045 

x,  metastable  lifetime 

10  ms 

<TPjpump  absorption  cross  section  (980  nm) 

a,  absorption  cross  section  (1.55  nm) 

T2,  polarization  relaxation  time 

60  fs 

Ld,  dispersion  length  (300  fs  pulse) 

450  cm 

Numerical  Methods 

The  nonlinear  Schrodinger  equation  was  integrated  using  the  split  step  Fourier 
transform  method.  This  is  an  example  of  a  spectral  algorithm.  These  are  described  by 
Strang23  and  DeVries.24  This  procedure  is  described  by  Agrawal4  and  a  FORTRAN 
program  for  the  algorithm  for  the  simplified  NLSE  without  higher  order  dispersion,  self- 
steepening,  or  Raman  scattering,  is  found  in  the  book  by  Newell  and  Moloney.25  This 
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program,  translated  into  the  C  programming  language,  is  the  basis  of  the  work  done  here. 
To  set  up  this  algorithm,  one  rearranges  Equation  2.35,  and  breaks  it  up  into  its  linear  and 
nonlinear  parts.  The  equation  is  then  expressed  as 


SU 


(2.39) 


where 


L  = 


a2  i  a3 

ST2  +  6  ST3 


(2.40) 


and 


NA  =  iy 


,2.  .  2i  a  a , ,2 , \  A  a|A| 


A  A  +  — — 
11  co05T 


(|a|2a)-tr 


A- 


3T 


(2.41) 


If  this  equation  is  formally  integrated,  it  yields 

•  U(4  +  h)  =  exp(|  £,)  exp|"  jN(z)dz  exp(^  £,]  U(^) , 


(2.42) 


V  5  J 

where  h  is  the  step  size.  The  operator  for  the  linear  step  is  carried  out  in  the  Fourier 
domain 


exp 


^hL(co)^ 
V  2  J 


B(^,t)  =  F” 


exp 


.  hL(co) 


(2.43) 


We  replace  the  operator  —  with  ico  in  Fourier  space  where  F  denotes  the  Fourier 

dl 

transform. 

The  nonlinear  operator  is  approximated  using  a  variation  on  the  trapezoidal  rule. 
For  the  first  iteration  the  operator  is  simply  assumed  to  be 

jN(z)dz  «  hN(z) . 

This  first  guess  is  used  in  Equation  2.42  to  derive  U(£+h).  This  new  value  of  U(^+h)  is 
used  to  derive  a  second  guess  for  the  nonlinear  operator.  Using  the  trapezoidal  algorithm 
for  integration  integration, 
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This  new  value  for  the  nonlinear  operator  can  be  used  to  form  a  new  guess  for  U(£+h). 
The  process  can  be  repeated  as  needed.  Four  or  five  iterations  are  sufficient,  though  for 
very  intense  pulses  as  many  as  10  were  used  for  the  results  reported  in  this  thesis. 

The  computer  program  developed  for  this  thesis  was  checked  in  a  number  of 
ways.  Initially  it  was  compared  to  published  results  by  other  groups.  This  process  is 
described  in  the  remainder  of  this  section.  In  addition,  some  results  were  compared  to 
those  of  another  student,  who  also  needed  to  simulate  a  similar  system,  and  who  had 
independently  written  a  computer  code  to  simulate  a  fiber  amplifier.  These  results  were 
found  to  be  identical.  In  general,  the  step  size  and  number  of  FFT  data  points  were 
increased  until  no  further  changes  were  observed.  For  very  intense  pulses,  as  many  as 
8 1 92  FFT  data  points  were  used,  though  often  1 024  points  proved  to  be  sufficient.  The 
width  of  the  time  window  of  the  simulation  was  varied  from  40  to  240  in  soliton  units, 
depending  on  the  application.  These  were  set  by  examining  the  results  in  time  and 
frequency  space  to  insure  that  the  intensity  of  the  pulses  had  fallen  to  very  small  amounts 
at  the  edges  of  the  time  and  frequency  windows.  For  a  pulse  with  an  amplitude  of  about  1 
in  soliton  units,  10  steps  per  dispersion  length  were  found  to  be  sufficient  for  accurate 
simulations.  As  many  as  10,000  steps  per  dispersion  length  were  used  when  simulating 
very  intense  pulses.  It  was  often  found  that  asymmetries  of  the  pulses  at  the  end  of  a 
simulation  indicated  a  problem.  This  was  inevitably  the  case  if  there  was  no  physical 
process  in  the  simulation  that  would  produce  this.  Finally,  wherever  possible,  analytical 
predictions  of  the  model’s  results  were  made.  Many  of  these  are  reported  in  the  thesis. 
Rough  agreement  between  the  numerical  and  analytical  results  also  encourages 
confidence  in  the  simulation.  Finally,  the  simulation  reproduces  many  phenomena 
observed  in  experiments.  These  will  be  pointed  out  in  latter  sections. 


Simulation  Results 


The  algorithm  described  above  was  carried  out  for  different  values  of  the 
constants  and  compared  with  published  results  to  verify  the  accuracy  of  the  program. 
These  results  also  illustrate  how  the  various  higher  order  effects  alter  the  soliton.  Figure 
2.1  was  the  result  of  propagating  a  fundamental  soliton  down  a  fiber  when  no  higher 
order  terms  are  present.  The  result  is  that  the  pulse  retains  its  hyperbolic  secant  shape 
indefinitely.  This  case  ignores  the  effect  of  intensity  loss,  but  otherwise  this  case  is 
approximately  correct  for  pulses  longer  than  about  one  picosecond. 

Figure  2.2  shows  the  effect  of  stimulated  Raman  scattering  (SRS).  The  case  given 

is  for  tr  =0.05  which  allows  for  direct  comparison  with  Figure  9.14  in  Agrawal.1  The 

parameter  Tr  can  range  from  6-8  fs.  The  effect  becomes  significant  when  pulse  duration 

approaches  this  value.  The  effect  of  SRS  is  a  phenomenon  known  as  self-frequency  shift, 

which  reduces  the  carrier  frequency  of  the  soliton.  This  causes  the  group  velocity  of  the 

pulse  to  be  reduced  in  the  anomalous  dispersion  regime.  Figure  2.3  shows  the  result  of 

the  phenomenon  known  as  self-steepening.  Here,  s=0.2,  which  makes  Figure  2.3 

comparable  with  Figure  5. 16  in  Reference  1 .  Again,  the  primary  effect  is  to  slow  down 

the  propagation  of  the  soliton  in  the  fiber.  Figure  2.4  shows  the  effect  of  higher  order 

dispersion.  Here  the  pulse  is  propagated  for  6  soliton  periods,  with  5=0.33  which  allows 
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for  comparison  with  the  paper  by  Liu  and  Puri. 
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Normalized  Intensity 


0  Dispersion  Lengths - 3  Dispersion  Lengths . 5  Dispersion  Lengths 


0  Dispersion  Lengths - 3  Dispersion  Lengths . 5  Dispersion  Lengths 


Figure  2.1  Soliton  pulse  shape  with  xR=0.0,  s=0.0,  5=0.0,  after  traveling  0,  3,  and  5 
propagation  distances.  Note  that  there  is  a  complete  overlap  of  the  lines 
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Figure  2.2  Pulse  shapes  after  0,  3,  and  5  dispersion  lengths.  tr=0.05.  The  Figure 
shows  how  the  pulse  slows  down,  due  to  Raman  shift. 
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Figure  2.3  Pulse  shapes  after  0,  5,  and  10  dispersion  lengths.  s=0.2.  Figure  shows 
how  the  pulse  slows  down  due  to  self-steepening. 
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Normalized  Intensity  Normalized  Intensity 


0  Soliton  Periods - 3  Soliton  Periods 


0  Soliton  Periods - 3  Soliton  Periods 


Figure  2.4  Pulse  shapes  after  0,  and  3  soliton  periods.  8=0.033.  Figure  shows  the 

effect  of  higher  order  dispersion. 
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Chapter  Three  :  Propagation  in  Birefringent  Optical  Fibers 

Normal  fibers  possess  random  linear  birefringence  due  to  strain  or  random 
imperfections  in  their  cores.  Linear  birefringence  of  optical  fibers  has  been  reviewed  by 
Kaminow.  The  effect  of  linear  birefringence  is  that  the  components  of  light  along  the 
axes  of  the  fiber  will  have  different  group  velocities.  This  tends  to  break  up  pulses.  In 
the  linear  regime,  power  oscillates  periodically  between  the  two  polarizations  with  a  beat 
length  Lb  ,  where  Lb  =  A./ An  and  An  is  the  difference  in  index  of  refraction  between  the 
two  axes. 

In  addition  to  the  linear  effect,  there  is  a  nonlinear  contribution  to  birefringence 
from  self  phase  modulation.28,29  The  polarization  dependent  nonlinearities  result  from  a 
modification  to  the  optical  Kerr  effect.  This  will  lead  to  a  change  in  the  nonlinear 
polarization  and  thus  to  an  altered  form  of  the  Nonlinear  Schrodinger  Equation.  The 
nonlinear  polarization  will  take  the  form  of 

Pnl(?,  t)  =  5  (xPx  +  yPy )  exp(-ico0t)  +  c.  c.,  (3.1) 


where  Px  and  Py  represent  envelope  functions  and  c.c.  denotes  the  complex  conjugate  of 
the  other  term  on  the  right  hand  side.  The  components  Px  and  Py  are  given  by 
3sf 


R  = 


■L,(x®)TE,EJE;+X««yE1E,E;+x<iEJEJE;), 


(3.2) 


where  the  indices  i  and  j  can  be  either  x  or  y.  The  functions  Ex  and  Ey  are  the  envelope 
functions  as  used  in  Equation  2.1.  Due  to  symmetry  relationships,  the  components  of  %(3) 
are  related  by 

3y(3)  (3)  (3)  (3) 

A  xxxx  A  xxyy  A  xyxy  A  xyyx  * 

If  we  use  this  relationship  in  Equation  3.2  and  3.1  for  the  polarization  at  frequency  co0  we 
get  the  x  component 
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,i(3E,E,E;+E,E,E;+EIExE;+EJE,E:)  = 


10.  v  (3) 

4  A/  xxxx 


eJ2+-|ev|2V+-(e*xev)ev 

xi  3 1  y I  )  *  3V  x  y;  y 


(3.3a) 


and  y  component 


(iv)"™" = ^xSl ± (JE^EyE; + e^e.e; +Eie,e; + eae;) = 

\  1  (3.3b) 

^x2{yE,f  +||E,rjEy  +i(E;E<)E>  . 

The  last  result  needs  to  be  related  to  an  effective  perturbation  of  the  refractive  index.  We 
can  relate  the  polarization  to  the  refractive  index  by 

Pi  =  £o£iEi  =  £o(n  +  Ani)2Ei>  (3-4) 

where  the  induee  i  indicates  either  the  x  or  y  axis.  The  last  expression  can  be 
approximated  by 

Pi  =  s0(n2  +2nAni)Ei.  (3.5) 


If  we  compare  3.5  to  3.3  we  see  that 


AnA=n„  l|Exf+^|E,|2k+i(E;Ey)E 


31  y|  J  x  3' 


(3.6a) 


AnyEy  =  nk|"f|Ey|2  +^|Ex|2lEy  +i(E',E  )E,  , 


where 


J  (3) 

n  = - y K  ’ 

k  o  A- xxxx  ‘ 

8n 


(3.6b) 


Equations  3.6  can  be  substituted  into  Equation  2.20  to  study  the  effect  of  polarization  on 
the  NLSE.  Equation  2.20  showed  that 
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k0  J  |F*AnF0dxdy 


•  W  WJ 

j  |F0*F0dxdy 

—00—00 

If  we  substitute  Equation  3.6  into  Equation  3.8,  we  get 

j  jF0*F0dxdy 

-00-00 

■  ]  J1F0 12 |ES |2 dxdy  +  ~  J  J1f0 I2 |Ej p dxdy  A,  +  ^]]|F,|2(i;Ej)Ejdxdy 


where  j  is  the  symbol  for  the  axis  orthogonal  to  the  i  axis.  If  we  expand  the  expressions 
for  the  electric  fields  we  find 


Ap,Ai  = 


00  00 

konk  j  j|F0|4dxdy 


F0*F0dxdy 


(3.10) 


|AiT  +||Ajf  Ai  +^(A*AJ)AJexp[-2i(pii  -p, 


If  we  group  variables  together,  we  get 


AP‘Ai  =^L{lAi|2  +f  [Aj|2  +^(AiAJ)Ajexp[-2i(Pii  “Pi 


(3.11) 


where  Aeff  was  defined  in  Equation  2.29.  The  last  equation  can  be  put  into  Equation  2.26 


which  was 


5A_  i  ^A_P1^A 

at  2  at2  “  6  at3  + 11 


(3.12) 


to  get 


^±  =  _r  ^_iR  a2A;  P3  a3Aj  q  A 
az  Pl  a  2P2  a2  6  a3  2Af 

.  •  I  A  |2  *  .  2  .  I  .  |2  1  .  /  .  *  .  \  . 


'iy\AifAi  +-iY|Ajf  Ai  +-iY(A*Aj)AJexp[-2i(pii  -pu)  . 


(3.13) 


The  Raman  effect30,  and  self-steepening  are  both  altered  by  birefringence.  Both 
of  these  terms  can  be  added  to  Equation  3.13.  The  Raman  term  will  be  examined  first.  It 
is  found  that  the  Raman  effect  has  the  general  form 

.  Su^t)  =u(t)£f](t-t')|u(t')|2dt'  + 

Raman 

(3.i4) 

v(t)|  f3(t-t')u(t')v*(t')dt'  + 

v(t)  |  f3(t - t')u * (t')v(t')exp(iR5z)dt , 

where  the  functions  u  and  v  are  the  normalized  envelopes  introduced  in  Equation  2.35.  In 
an  isotropic  medium,  fi=f2+  f3+  U-  It  is  also  found,  in  general,  that  f2=  f3=  U-  This  means 
f 

that  f2  =  — .  The  effect  of  the  last  term  can  be  ignored,  because  the  phase  mismatch 

between  the  terms  will  make  them  average  to  zero.  For  pulses  longer  than  about  lOOfs 
we  can  approximate  the  response  functions  by  the  first  term  of  their  Taylor  expansions  as 
was  done  in  the  derivation  of  Equation  2.32.  This  leaves  us  with 

,5u^t)  =Tr^u(t)f+Iu(t)|v(t)|2).  (3.15) 

Raman  V  J  ' 

The  Raman  and  self-steepening  terms  can  be  added  to  Equation  3.13.  The  last 
term  of  3.13  can  be  neglected  as  it  tends  to  be  rapidly  varying.  The  soliton  unit 
normalization  introduced  for  Equation  2.36  can  also  be  used  here.  The  result  is 


.  5u£t)  =T( u(t)|u(t)|2  + 1 U(t)| v(t 


(3.15) 


.  3u  Su  j  _ 

1  52u  f,  ,2  2,  ,2^ 

- T+  U  +— V  u  + 

2  dx 2  \n  3  / 


(3.16a) 


0|uf 

1  3v 2 

| 

dx 

3  dx 

U-1S -  U  UH - V  U 
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.  dv  dv 

-l  —  +  6  — 

Kdt,  d % 


1  S2v  ( I  ,2  2,  ,2^1 

- r-+  V  +— U  V  + 


T  [M.  +  I® 

Ri  dx  3  dx 


(3.16b) 


R  +  v-is—  v  v  +  -  u  v  . 


7i  c  An 

Here  8  =  ,  where  An  is  the  difference  in  index  of  refraction  between  the  two 

D(a)a,0 

axes.  When  higher  order  terms  are  neglected,  this  reduces  to 


du  d  u 


d  v  dv 
_1  +  ^  = 
\dZ,  d^J 


1  d2  u 

2  dx  2 


2  +  M2  +  “|v|2  u 


(3.17) 


15  v  ,  ||  I2  ,  2  I  I2  ] 

+ilul  Jv- 


For  a  300fs  soliton,  normal  values  of  5  will  range  from  0.02  to  0.20.  Fibers  with 
increased  birefingence  may  be  used  in  order  to  eliminate  the  effect  of  random  fluctuations 
in  the  birefringence. 

The  question  now  arises  as  to  whether  Equation  3.17  has  soliton  solutions,  and  if 
it  does  how  they  behave.  In  fact,  a  wide  variety  of  soliton  solutions  have  been  found. 

This  question  has  been  examined  by  Trantnik  and  Sipe31,  Christodoulides  and  Joseph32, 
Tasgal  and  Potasek33,  and  by  Haelterman  et  al.34  One  form  which  has  been  quite  useful 
was  found  by  Ueda  and  Kath.35  For  the  case  of  equal  amplitudes  on  both  birefringent 
axes,  they  found  the  solution 


30 


(3.18) 


u  =  r)  sec  h 


2 

1+V 


1/2 


exPi  i 


2)r\2  , 

1  +  —  —  +  —  E  -  i5x 
3/  2  2  b 


v  =  r|  sec  h|  1 1  + 


\l/2 


3) 


px 


exp  I  i 


v1+  3 ) 


T1  O  , 

y  +  Y^  +  iSx, 


The  Equations  3.18  state  that  the  vector  soliton  should  be  narrower  by  (5/3)1/2.  The  phase 
shift  as  a  function  of  distance  will  also  be  altered.  Recently  a  more  general  approximate 
soliton  form  was  found  by  Afanasjev.36  The  form  can  have  an  arbitrary  polarization 
angle  and  allows  study  over  a  wider  range  of  pulse  input  angles.  Since  the  results  have 
only  recently  become  available,  they  were  not  used  in  this  project. 

The  nonlinear  term  has  two  primary  effects.  One  is  referred  to  as  soliton 
dragging.37’38  If  the  pulses  have  an  intensity  over  a  certain  limit39,  the  light  on  the  fast 
and  slow  axes  is  shifted  in  frequency  so  that  they  travel  at  the  same  group  velocity.  This 
effect  is  shown  in  Figure  3.1.  It  shows  the  result  of  the  simulation  of  propagating  a 
soliton  of  unit  normalized  intensity  5  soliton  periods  in  a  birefringent  fiber  without  higher 
order  terms.  In  this  figure  the  input  amplitude  in  both  polarizations  is  equal.  This  effect 
can  eliminate  the  problem  of  pulse  break  up  due  to  linear  birefringence.  If  the  pulses  do 
not  have  equal  amplitudes  along  both  axes  the  behavior  changes.40  If  the  pulses  are 
intense  enough  that  the  dispersion  length  is  less  than  a  beat  length  the  pulse  is  stable, 
though  energy  can  be  exchanged  between  the  polarizations.  If  the  beat  length  is  less,  the 
energy  on  the  fast  axis  of  the  fiber  is  transferred  to  the  slow  axis. 

The  frequency  shift  produced  by  soliton  dragging  can  be  used  to  make  a  switch.11 
The  switch  works  in  the  following  fashion.  If  a  pulse  is  launched  into  a  piece  of 
birefringent  fiber  along  one  of  the  birefringent  axis,  it  emerges  from  the  other  end  with  its 
frequency  unaltered.  If,  however,  a  second  pulse  on  the  other  axis  is  launched  along  with 
it,  the  frequency  of  the  pulses  will  be  altered.  A  filter  is  placed  at  the  output  such  that  it 
will  not  transmit  the  unshifted  light  but  will  transmit  the  pulse  if  it  is  shifted  by  soliton 
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dragging.  This  combination  will  act  as  an  And-gate  since  it  will  only  transmit  a  pulse  if 
both  pulses  are  present. 

A  second  useful  manifestation  of  nonlinear  birefringence  is  intensity-dependent 
polarization  rotation.41  The  polarization  of  light  which  is  originally  polarized  along  one 
of  the  birefringent  axes  will  not  change.  This  can  be  seen  in  Equation  3.16.  The 
polarization  will  also  not  change  if  the  pulse  originally  has  equal  components  along  both 
axes.  In  this  case,  the  nonlinear  change  in  the  index  of  refraction  is  the  same,  so  it  does 
not  alter  the  beat  length.  It  should  be  noted,  however,  that  the  factor  of  2/3  in  the  Ken- 
term  means  that  the  self-phase  modulation  for  a  pulse  with  equal  components  along  the 
birefringent  axes  will  be  less  than  that  for  a  pulse  which  is  polarized  along  one  axis.  This 
will  have  significance  in  the  study  of  birefringent  fiber  devices.  For  other  polarizations, 
the  polarization  of  light  will  evolve  in  a  fashion  which  depends  upon  its  intensity.  This 
effect  is  greatest  if  the  light  is  polarized  at  an  angle  re  /  8  from  one  of  the  birefringent 
axes.  Intensity-dependent  polarization  rotation  has  been  used  to  construct  switches  and  to 
mode-lock  lasers. 
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Normalized  Intensity  Normalized  Intensity 


Slow  Axis 


Fast  Axis 


Slow  Axis  - Fast  Axis 


Figure  3.1  Shows  the  effect  of  pulse  propagation  in  5  soliton  periods  in  a 
birefringent  fiber.  The  birefringence  is  5=0.5,  The  input  pulse  is  a  unit  amplitude 
vector  soliton  with  equal  amplitude  in  both  polarizations. 
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Chapter  Four :  Pulse  amplification 

Having  studied  the  propagation  of  pulses  in  a  fiber,  we  will  now  examine  the 
propagation  of  pulses  in  fibers  which  act  as  a  gain  medium.  The  fiber  amplifier  is  one  of 
the  major  devices  in  a  fiber  laser.  This  chapter  describes  previous  work  on  this  topic. 

The  basic  equation  which  I  used  to  simulate  the  fiber  amplifier  is  derived.  While  this 
equation  is  not  new,  I  have  not  been  able  to  find  a  satisfactory  derivation  in  another 
source.  Having  done  this,  I  report  on  novel  work  I  have  performed  which  describes  the 
interaction  of  the  effect  of  birefingence  and  pulse  break-up  due  to  the  amplifier. 


Pulse  Simulation  Theory 


*  The  basic  properties  of  erbium  doped  fibers  have  been  reviewed  by  Ainslie42. 

Early  experimental  results  were  reported  by  Nakazawa  et  al.43,  Ainslie  et  al.44,  and 
Khrushchev  et  al.45  General  descriptions  of  fiber  amplifiers  can  be  found  in  Giles  et  al.46 
and  Desurvire.47  There  are  two  lines  of  investigation  that  have  been  pursued,  to  study  the 
amplification  of  pulses.  Both  start  by  using  Equation  2.35  to  describe  how  the  electric 
field  evolves  in  the  fiber.  In  addition  to  this,  one  needs  equations  that  describe 
amplification.  Mel'nikov  et  al  48,  Wang  and  Chi49,  Maimistov50,  Gabitov51,  Nakazawa52, 
and  Gross  and  Manassah53  start  with  the  Bloch  equations3, 


da 


21 


dt 


=  i(co  -co0)a21  +^L(Pll  -p22)-^r- 


2  h 


d  („  „  (Pii  P22)  (pn  P22 )o 

dt(Pn  P22)-  (a2i  CT2lj 


(4.1) 


(4.2) 


In  these  equations  p  is  the  dipole  moment  in  the  interaction  Hamiltonian,  EQ  is  the  field 
amplitude  of  the  illuminating  light,  and  (oQ  is  the  resonance  frequency  for  the  transition. 
The  variable  T2  is  the  decay  time  of  the  polarization  and  x  is  the  population  decay  rate. 
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The  variables  p  are  the  elements  of  the  density  of  states  matrix. 

slowly  varying  part  of  p,  which  is  related  to  G21  by 

The  variable  (J21  is  the 

P21W  =  cr21(t)exp(-icot) 

(4.3) 

where 

CT21  =  ai2  • 

(4.4) 

The  mean  dipole  moment  of  the  atom  is 

^p(t))  =  2p[Rea21(t)cos(o)t)  +  Ima21(t)sin(cot)] . 

(4.5) 

We  can  relate  this  to  the  mean  polarization  in  the  material  by 

P  =  N(n), 

(4.6) 

where  N  is  the  density  of  atoms.  The  population  inversion  per  unit  volume  is  given  by 
AN  =  N(pu  -  p22)  •  (4-7) 

If  we  put  Equations  4. 3-4.7  into  Equations  4.1  and  4.2,  we  see  that  they  form  together 
with  Equation  2.35  a  set  of  three  Equations  with  the  variables  E,  P  and  AN.  This  set  of 
equations  can  be  set  up  as  an  initial  value  problem  and  solved  to  model  the  amplification 
of  a  pulse.  The  Bloch  equation  approach  accurately  captures  the  physics  of  amplification 
and  can  be  used  to  study  very  short  pulses.  The  problem  is  that  this  method  takes  a  great 
deal  of  computation  and  so  another  method  is  used  where  it  is  possible  to  do  so. 

This  simpler  approach  has  been  used  by  Hodel  et  al.54-55,  Gagnon56,  Belanger57, 

Nakazawa58,  Romagnoli59’60,  and  Agrawal61>62>63.  If  the  pulses  are  longer  than  the 
parameter  T2  which  is  on  the  order  of  100  fs;  then  this  simpler  approach  can  be  used. 
This  is  the  case  for  all  but  the  shortest  pulses.  This  method  has  been  used  for  the 
simulations  described  in  this  thesis.  Here  the  gain  is  approximated  by 


g(®)  =  gp[1-(®  -coo)  T2 
where  gp  is  the  peak  gain.  This  is  sometimes  referred  to  as  the  parabolic  gain 
approximation.  The  time  dependence  of  the  peak  gain  is  given  by 


(4.8) 
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(4.9) 


3gp_g0-gp  gP|A| 


These  equations  can  be  derived  from  the  Bloch  equations,  when  the  derivatives  are  not 
too  large.  One  starts  by  assuming  that  the  derivative  in  Equation  4.1  is  zero.  Equation 
4.1  and  its  complex  conjugate  are  linearly  independent.  The  sum  and  difference  of 
Equation  4.1  and  its.  complex  conjugate  are  taken.  This  leaves  us  with 


(a2i  -CT*2i)  =  iT2(®-®o)(^2.  +  °*21) 


(4.10a) 


(a21  +  a  *21 )  =  iT2  (co  -  03  0)(a2l  -  a  *21 )  +  -^E°  (Pn  -  p22 ) . 

When  the  term  involving  the  sum  of  the  cr2i  terms  is  eliminated  one  is  left  with 
(  \  _  ^pEp  (Pll~p22) 

^  21  21 ’  ft  1  +  T22(co  -  (o0)2 

Substituting  this  into  Equation  4.2,  we  get 

—  (  _  )  _P2T2E02  (pn~P22)  _  (Pn  ~P22)~(Pn  ~P22)0 

dtlPn-P22j  =  -  h2  1  +  J2^_moy- 


(4.10b) 


(4.11) 


(4.12) 


The  normalized  line  shape  function  in  the  dipole  moment  model  is  given  by  Yariv9  in 
Equation  8.1-20  as 


g„(®)  = 


1  +  (co-co0)2T22 


(4.13) 


Substituting  this  into  Equation  4.12,  one  obtains 
d  /  \  P  gn(®)E0  (pll  “P22)  (Pll  _  P22 )  (pl  1  “  P22)0 

dt  P’l  P22/  2h2  ~  "  • 

The  gain  parameter  is  defined  as 

g(co)  =  AN  a(co)  =  AN  agn(ffl)  =  gpgn(co) , 

where  a  is  the  absorption  cross-section.  Clearly,  the  prefactor  in  Equation  4.8  is 


(4.14) 


(4.15) 
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(4.16) 


gP  =  AN  a. 

In  Reference  9,  Equation  8.2-7  states  the  gain  parameter  is  expressed  as 

,  X  kx"(©) 

g((Oj  =  - - 2 - » 

n 


(4.17) 


where  y”  is  the  imaginary  part  of  the  atomic  susceptibility.  It  can  be  expressed9  as 


X"M  =  :TTANgnM- 


(4.18) 


When  the  last  two  equations  are  combined  we  get  an  expression  for  peak  gain.  When  this 

is  compared  to  Equation  4.15  one  sees  that 

cop2  cop2n 
2cmnh  2c  zh .  ’ 


where  the  fact  that  n2  =  —  was  also  used.  When  4.14  is  multiplied  by  Nor,  the  result  is 

so 

dgP  _  E2gn(m)Eo2gP  gp-go 
dt  2  fi2  x 

The  last  equation  can  be  expressed  in  terms  of  the  variable  A,  which  is  described  by  the 
square  root  of  the  power,  instead  of  E0.  This  can  be  done  using  the  relationship 


C£  It?  I2 

- E0 

2nA 


=,=W1. 


We  are  also  interested  in  using  the  time  average  of  the  electric 


field.  This  can  be  expressed  as  E  =  —  |E0|  .  When  these  expressions  are  substituted  into 


the  last  equation  the  result  is 

dgp  _  E2  1  Ap  ngp  gp  -g0 
dt  2ft2Aeffcs  x 

When  a  factor  of  a  is  factored  out  of  the  last  expression,  the  result  is 

dgp  |Al2  gp  _gp  -g0 

dt  Es  t 
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Which  is  the  expression  in  Equation  4.9.  If  x  is  much  greater  than  the  pulse  width,  we 
can  ignore  the  last  term  on  the  left  hand  side.  If  we  integrate  we  get 


gP(T)  =  go  expf~S  £jA|2dx  . 


(4.19) 


This  last  equation  allows  us  to  include  the  effects  of  saturation  which  is  expressed  by  the 
parameter  S.  Here  the  saturation  energy  Es  is 
fi  conA„ 


£  _  xeff 

G 


(4.20) 


where  a  is  the  transition  cross  section.  If  we  change  (co-co0)  to  the  operator  i—  in  the 

dt 

time  domain,  we  get  that  the  form  for  a  in  Equation  2.35  becomes 


a  =  gP(t) 
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1  +  T. 


2  ax2 


dt1 


(4.21) 


In  normalized  units,  this  becomes 


an  =g0LDexp[-s£|U|2dx 


f  2  d 2  " 
2  dx\ 


1  +  T 


(4.22) 


where  x2  T2/TQ  and  S  Po  TQ/Es  The  saturation  energy  is  about  10  pJ.  The  parameter  S 
should  then  be  around  10‘6. 

Since  amplification  in  an  erbium  amplifier  is  largely  independent  of  polarization, 
the  extension  of  Equation  4.13  to  the  case  of  two  polarizations  is  straightforward.  One 
must  remember  that  both  polarizations  are  contributing  to  polarization,  leaving  an 
amplification  of  the  form 


«n  =goLDeXP 


-s£K+M>]v 


1  +  x 


2  & 


,2  X 


dx2J 


(4.23) 


Simulations  of  this  kind  have  been  carried  out  by  Yang  and  Wang64.  They  included  the 
effects  of  Raman  response  and  self-steepening.  They  demonstrated  that  light  along  the 


two  axes  tend  to  couple  and  that  the  spacing  between  them  is  reduced  as  the  pulses  are 
amplified. 

The  equations  shown  here  treat  only  the  longitudinal  effects  in  a  fiber  amplifier. 
This  is  justified  in  reference  4.  It  has  been  found  that,  in  general  combining  transverse 
with  longitudinal  effects  makes  only  a  small  difference.  The  importance  of  studying 
transverse  mode  and  doping  shapes  is  of  importance  to  optimizing  fiber  amplifiers.  This 
has  been  done  by  Desurvire  et  al.  ,  Ohashi  ,  Pedersen  et  al.  ,  and  Pfeiffer  and  Bulow  . 
The  results  reported  show  the  importance  on  concentrating  the  light  from  the  pump,  the 
signal  to  be  amplified,  and  the  dopant  ions  in  the  center  of  the  core.  For  this  reason,  fiber 
amplifiers  often  have  smaller  cores  than  normal  telecommunications  fibers. 


Simulation  Results 

I  simulated  the  propagation  of  fundamental  solitons  with  a  pulse  width  of  300 
femtoseconds  introduced  into  fiber  amplifiers  with  a  gain  of  10  dB  per  dispersion  length. 
Here  t2=0.2.  Numerically,  these  were  the  most  difficult  simulations  performed  due  to 
the  very  high  intensities.  Results  for  the  amplifier  without  birefringence  were  compared 
with  those  of  another  student  who  independently  developed  a  similar  program.  The 
results  were  the  same.  The  width  of  the  time  window  of  the  simulation  and  the  number 
of  FFT  points  were  increased  until  the  pulse  was  contained  in  the  window  in  both  time 
and  frequency  space.  The  time  window  was  240  soliton  widths  wide  and  8192  FFT 
sample  points  were  used.  Due  to  the  high  intensities  10,000  steps  per  dispersion  length 
were  taken  in  the  simulation.  Increasing  the  number  of  steps  did  not  alter  the  results  so 
they  were  assumed  to  have  stabilized. 

The  results  can  be  seen  in  Figures  4.1  -  4.2.  The  effect  of  gain  dispersion  is  to 
limit  the  height  of  the  pulse  and  to  cause  it  to  break  up.  The  effect  seen  in  both  time  and 
frequency  space  is  that  side  peaks  are  continually  built  up  after  the  pulse  reaches  a  critical 
energy.  As  the  peaks  gain  intensity  they  are  compressed  in  time.  This  causes  them  to 
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expand  in  bandwidth  and  so  eventually  two  smaller  pulses  will  be  energetically  favored 
over  one  large  one. 

The  effect  of  gain  saturation  can  be  seen  in  Figure  4.3.  Here,  everything  is  the 
same  as  with  the  last  Figures  except  that  the  saturation  constant  is  set  to  S=0.01 .  This 
was  done  to  illustrate  the  effect  of  saturation.  As  shown  above,  the  actual  saturation  is 
much  smaller.  The  saturation  tends  to  reduce  the  maximum  intensity  of  the  pulse,  and  to 
cause  the  leading  edge  to  be  built  up  relative  to  the  trailing  edge. 


3  Dispersion  Lengths - 5  Dispersion  Lengths 
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Normalized  Time 


Figure  4.1  Shows  how  a  pulse  is  modified  in  an  amplifier.  The  gain  is  10  dB  per 
dispersion  length.  There  is  no  saturation.  The  results  are  given  in  terms  of 

normalized  time 
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3  Dispersion  Lengths - 5  Dispersion  Lengths 


Normalized  Frequency 


Figure  4.2  Shows  the  how  a  pulse  is  modified  in  an  amplifier.  The  gain  is  10  dB  per 
dispersion  length.  There  is  no  saturation.  The  results  are  given  as  a  function  of 

frequency. 
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intensity 


Figure  4.3  Shows  the  effect  of  pulse  amplification.  The  gain  is  10  dB  per  dispersion 
length.  Results  are  given  both  with  no  saturation,  and  for  S=0.01.  The  amplifier  is 
3  dispersion  lengths  long.  The  results  are  given  in  terms  of  normalized  time 
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time 


. 1  dispersion  length  - 2  dispersion  lengths 

Figure  4.5  Result  of  pulse  propagation  in  an  amplifier.  The  gain  is  10  dB  per 
dispersion  length.  There  is  no  saturation,  and  5=0.0. 

The  result  of  polarization  was  examined  in  a  series  of  simulations.  The  amplifier 
parameters  were  the  same  as  in  the  last  paragraph.  Figures  4.4  and  4.5  show  the  case  in 
which  8=0.0,  where  8  is  the  parameter  measuring  birefringence  given  in  3.16.  In  this  case 
the  pulse  had  equal  intensities  on  both  axes  of  the  fiber.  The  results  for  both  axes  were 
the  same.  After  one  dispersion  length  the  pulse  has  been  amplified  and  compressed.  The 
phase  is  nearly  constant  over  the  high  intensity  part  of  the  pulse,  which  is  expected  for  a 
soliton.  The  result  is  that  after  propagating  two  dispersion  distances  the  pulse  begins  to 
break  up  as  in  the  case  in  which  we  neglected  polarization.  Before  pulse  breakup  the 
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phase  has  a  single  maximum,  but  after  breakup  there  is  more  modulation  and  multiple 
maxima.  The  maxima  and  the  central  flat  region  of  the  phase  correspond  to  the  location 
of  the  pulses.  Again  they  are  acting  like  solitons  in  that  they  have  a  uniform  phase. 

The  effect  of  introducing  birefringence  is  shown  in  Figures  4.6-4. 8.  The 
birefringence  was  altered  so  that  8=0.5.  The  input  pulse  had  an  angle  of  7i/8  to  the  slow 
axis  of  the  fiber.  The  light  in  the  two  polarizations  is  beginning  to  separate,  most  notably 
in  that  light  in  the  fast  axis  is  beginning  to  build  up  in  a  broad  peak  in  the  front  of  the 
pulse.  Clearly  the  polarization  varies  within  the  pulse.  The  phase  difference  between  the 
two  polarizations  is  beginning  to  increase  and  it  varies  by  a  great  deal  across  the  pulse  as 
well.  Figure  4.8  can  be  thought  of  as  the  effect  of  placing  a  polarizer  at  the  output  of  the 
fiber.  The  amplifier  and  polarizer  working  together  have  the  effect  of  introducing 
significant  distortions.  Figure  4.9  shows  the  effect  of  introducing  saturation  into  the 
amplifier.  Here,  the  saturation  is  S=0.01.  The  pulse  has  experienced  less  amplification 
though  this  effect  is  much  less  pronounced  for  the  light  on  the  fast  axis. 
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slow  axis - fast  axis 


Figure  4.6  Result  of  pulse  propagation  in  an  amplifier.  The  angle  of  the  input  pulse 
to  the  slow  axis  of  the  fiber  is  71/8.  The  gain  is  10  dB  per  dispersion  length.  There  is 

no  saturation,  and  5=0.5. 
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Figure  4.7  Result  of  pulse  propagation  in  an  amplifier.  The  angle  of  the  input  pulse 
to  the  slow  axis  of  the  fiber  is  71/8.  The  gain  is  10  dB  per  dispersion  length.  There  is 

no  saturation,  and  8=0.5. 
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Figure  4.8  Components  of  the  pulse  shown  in  Figure  4.5  either  in  the  input 
polarization,  or  perpendicular  to  the  input  polarization,  as  shown. 
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- Original  Polarization . Orthoginal  Polarization 

Figure  4.9  Result  of  pulse  propagation  in  an  amplifier.  The  angle  of  the  input  pulse 
to  the  slow  axis  of  the  fiber  is  71/8.  Pulse  amplification  gain  is  10  dB  per  dispersion 
length.  The  saturation  is  S=0.01.  The  birefringence  is  5=0.5. 

Fiber  amplifiers  are  a  key  component  that  have  made  possible  a  number  of  major 
advances  in  optical  systems.  In  particular,  they  have  made  fiber  lasers  practical.  A  good 
model  of  pulse  amplification  in  fiber  amplifiers  was  needed  to  build  a  good  model  of  the 
figure  eight  laser.  A  model  was  described  here.  It  was  first  applied  to  the  case  of  the 
fiber  laser  without  birefringence.  This  work  was  verified  by  comparing  it  to  the  results  of 
another  student  who  was  independently  developing  a  model  of  the  fiber  amplifier.  I 
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carried  out  novel  research  in  the  area  of  pulse  amplification  in  the  birefringent  amplifier. 

I  showed  that  pulses  tend  to  break  up  in  these  amplifiers,  as  they  do  without 
birefringence.  It  is  also  known,  as  mentioned  in  Chapter  3,  that  the  components  of  a 
pulse  which  are  along  different  birefringent  axes  of  the  fiber,  tend  to  couple  to  each  other. 
I  further  showed  that  this  coupling  phenomenon  and  the  pulse  break-up  effect  coexisted. 
This  meant  that  the  amplifier  produced  a  train  of  coupled  fibers.  This  effect  may  explain 
the  fact  that  fiber  lasers  have  been  observed  to  produce  bursts  of  pulses  when  their 
amplifiers  are  being  pumped  intensely.  In  carrying  out  these  investigations,  the 
simulation  for  this  essential  component  of  a  laser  was  thoroughly  verified. 
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Chapter  Five  :  Nonlinear-Optical  Fiber  Loop  Switches 


The  component  of  a  figure  eight  laser  which  performs  passive  mode-locking  is  the 
Nonlinear  Optical  Loop  Mirror  (NOLM).  Because  of  its  critical  role  in  the  laser,  I  carried 
out  very  extensive  simulations  of  this  device.  After  duplicating  the  previous  reported 
work  of  others,  I  carried  out  novel  work  in  simulating  how  optical  pulses  were  effected  by 
an  NOLM  constructed  of  birefringent  fiber.  The  aspects  of  the  loop  studied  were  those 
that  effect  the  operation  of  figure  eight  lasers.  In  order  to  be  able  to  predict  the  pulse 
energies  and  widths  produced  by  a  laser,  one  needs  to  know  the  energies  for  which  pulses 
are  transmitted  by  the  NOLM.  Since  the  behavior  of  the  laser  is  sensitive  to  the  pulses 
being  distorted,  the  effect  of  the  loop  on  the  pulse  width  was  considered.  The  low 
intensity  transmission  of  the  loop  was  examined  because  this  will  be  the  state  of  the 
system  when  the  laser  is  starting  up.  One  must  also  know  how  the  operation  of  the  loop 
mirror  will  be  effected  as  the  parameters  of  its  components  are  altered.  Specifically, 
experimentally  a  great  deal  of  effort  is  always  spent  adjusting  polarization  controllers  to 
optimize  the  behavior  of  the  laser,  so  I  have  examined  the  impact  of  the  polarization 
controller  on  the  behavior  of  the  switch.  I  was  able  to  develop  a  simple  analytical  model 
which  allows  for  the  approximate  prediction  of  pulse  energy  required  for  switching  in 
these  devices.  Simple  linear  analysis  of  the  system  proved  to  be  enough  under  some 
conditions  to  predict  the  low  intensity  transmission.  These  models  shed  light  on  the 
physics  of  the  systems  and  give  the  ability  to  do  quick  analysis  for  laser  design. 


Previous  Work 

I  will  summarize  previous  work  on  these  devices  and  I  will  mention  at  least  some 
of  their  applications.  Nonlinear  fiber  switches  were  first  described  by  Doran  and  Wood.69 
The  basic  idea  for  a  nonlinear  switch  of  this  type  was  proposed  by  Otsuka.9  The 
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Nonlinear-Optical  loop  mirror  works  by  using  the  principle  of  a  nonlinear  Sagnac 
interferometer,  which  is  shown  in  Figure  5.1 .  Light  is  passed  through  one  input  port  of  a 
directional  coupler  which  acts  as  a  beam  splitter.  The  splitting  ratio  is  chosen  to  be 
something  other  than  50/50.  The  output  ports  of  the  directional  coupler  are  then  formed 
into  a  loop  of  fiber.  This  means  that  the  input  light  is  formed  into  two  counter- 
propagating  beams  in  the  loop.  Due  to  the  nonlinearity  of  the  fiber,  one  beam  will 
experience  a  greater  phase  shift  than  the  other,  since  the  beams  do  not  have  the  same 
amplitude.  The  intensity  difference  will  be  proportional  to  the  input  intensity.  This  phase 
difference  causes  switching  when  the  two  pulses  recombine  at  the  directional  coupler.  If 
the  beams  are  in  phase,  they  are  passed  out  the  port  they  came  in.  Effectively  they  are 
reflected.  If  they  are  n  radians  out  of  phase,  the  light  is  switched  to  the  other  port,  so  they 
are  transmitted.  The  sensitivity  of  the  device  can  be  greatly  increased  by  inserting  a  fiber 
Bragg  grating  into  it.70  This  may  be  useful  for  demultiplexing  applications. 

A  second  form  of  this  device  is  the  Nonlinear  Amplifying  Loop  Mirror 
(NALM).  This  device  employs  a  different  means  of  producing  beams  of  different 
intensities.  In  this  device  a  fiber  amplifier  is  placed  at  one  end  of  the  loop.  The 
directional  coupler  can  have  a  50/50  splitting  ratio.  If  the  amplifier  is  not  saturated,  light 
traveling  in  one  direction  is  amplified  at  the  beginning  of  its  passage  through  the  loop  and 
experiences  a  greater  phase  shift  than  light  going  in  the  other  direction  since  this  light  is 
amplified  only  after  traversing  the  loop.  The  phase  difference  is  caused  by  the  fact  that 
the  pulses  pass  through  most  of  the  fiber  in  the  loop  with  different  pulse  energies.  This 
phase  difference  again  produces  intensity  dependent  switching. 


s 
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Figure  5.1  Nonlinear  Optical  Loop  Mirror 


These  devices  have  been  used  for  a  number  of  applications.  The  switching  of 

79 

solitons  in  them  was  demonstrated  experimentally  by  Blow,  Doran,  and  Wood.  The 
NOLM  can  be  used  for  many  applications  in  optical  communications.  Their  use  in 
filtering  noise  in  optical  systems  has  also  been  demonstrated.74,75  When  cw  light  of  two 
different  wavelengths  is  introduced  into  such  a  switch,  multigigabit  pulse  trains  can  be 
produced.76  When  properly  configured  they  can  be  used  for  the  adiabatic  amplification  of 
solitons.77  It  has  also  been  studied  as  a  mode-locking  device  in  fiber  lasers.78 


Theory 

Most  previous  work  has  analyzed  these  devices  in  terms  of  their  cw  performance. 
This  work  will  be  briefly  described  here  since  it  helps  to  shed  light  on  how  these  devices 
work.  Take  Ei  to  be  the  cw  input  electric  field.  The  variable  a  is  the  power  transmission 
of  the  direction  coupler.  The  electric  fields  of  the  signals  in  the  loop  become 
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(5.1a) 


E3  =  a 1/2  Ei 
and 


E4  =  i(l-a)1/2Ei.  (5.1b) 

Where  E3  and  E4  are  the  two  counter  propagating  beams.  In  the  optical  fiber,  light  is 
subject  to  the  optical  Kerr  effect,  so  that  the  phase  is  shifted  by 


27ink|E|2L 

X 


(5.2) 


If  the  input  power  is  not  split  evenly  between  the  two  directions,  the  light  will  undergo  a 
different  phase  shift  as  it  travels  in  each  of  the  two  directions.  This  will  cause  intensity 
dependent  switching  when  the  light  is  recombined.  The  transmitted  light  will  be  given  by 


Eout  =  a1/2E3e^+i(l-a)1/2E4  =aE,  e*-(l-a)E,  = 
{ [1  +  cos(<(»)  +  i  sin(<j))]a  -  1}E,  , 


(5.3) 


where  <j)  is  a  relative  phase  shift  produce  in  E3.  After  the  light  passes  the  loop  the  output 
power  is  given  by  taking  the  magnitude  of  Equation  5.3  and  rearranging. 


E 

°l'*2  =  { [1  +  cos(<t>)  +  i  sin(<j>)]a  - 1}  { [1  +  cos(<j>)  -  i  sin((())]a  - 1}  = 

lEi| 

1  -  2a(l  -  a)[l  +  cos((j))] 


(5.4) 


If  the  relative  phase  shift  is  produced  by  the  Kerr  nonlinearity,  in  the  cw  case  we  get 


eout|2=|Ein|2  l-2a(l-a)<  1  +  cos 


(l-2a)|E,  \l2nn*L 
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(5.5) 
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where  L  is  the  length  of  the  loop.  By  examining  Equation  5.5,  the  general  features  of  the 
action  of  the  NOLM  can  be  found.  The  phase  shift  is  clearly  dependent  on  the  intensity 
of  the  input.  The  intensity  needed  for  switching  can  be  reduced  in  a  number  of  ways. 

The  loop  can  be  made  longer,  which  will  allow  more  interaction  path  over  which  the 
phase  shift  will  be  built  up.  The  splitting  ratio  of  the  directional  coupler  can  be  made 
larger.  This  will  mean  that  the  difference  in  nonlinear  phase  shift  per  unit  length  will  be 
larger.  Finally,  one  could  use  a  material  with  a  different  value  of  nk,  however,  this  is  not 
generally  practical.  These  general  features  carry  over  to  the  study  of  how  pulses  act  when 
introduced  into  an  NOLM. 

It  has  been  found  that  hyperbolic  secant  shaped  pulses,  those  with  an  intensity 
envelope  function  of  the  form 

Asech2(t), 

can  be  almost  completely  transmitted  through  an  NOLM  for  a  correct  value  of  the 
amplitude  A.  It  should  be  noted  that  unlike  the  situation  with  a  soliton,  the  width  and 
height  of  the  hyperbolic  secant  shaped  pulse  are  not  coupled.  Pulses  of  this  form  can  also 
form  solitons  if  they  have  enough  energy.  Solitons  have  the  property  that  they  have 
uniform  phase  across  their  temporal  extent.  After  the  input  pulse  is  split  by  the 
directional  coupler,  the  counter  propagating  pulses  will  form  themselves  into  solitons 
accompanied  by  low-intensity  pedestal.  The  pedestal  is  referred  to  as  dispersive  wave 
radiation.  The  rate  at  which  the  phase  changes  will  depend  on  the  energy  of  the  soliton, 
hence  a  phase  difference  builds  up  between  the  two  counter-propagating  pulses. 
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Simulation  Results:  Nonbirefringent  NOLM 

The  transmission  of  hyperbolic  secant  shaped  pulses  was  extensively  studied. 

The  result  of  simulations  of  NOLMs  in  which  the  effect  of  birefringence  is  not  considered 
will  be  described  here.  These  results  have  been  already  described  in  length  by  other 
researchers.  The  results  are  given  here  because  their  replication  gives  confidence  that 
my  simulations  are  working  well.  The  discussion  of  these  results  will  also  make  it  easier 
to  think  about  the  results  for  the  birefringent  NOLM.  In  addition,  the  model  for 
predicting  the  pulse  energy  needed  for  maximum  transmission  is  new.  Hyperbolic  secant 
shaped  pulses  of  varying  amplitude  were  introduced  into  a  simulated  NOLM.  First,  the 
pulses  were  split  into  two  subpulses  using  the  relationships  in  Equation  5.1.  The 
propagation  of  the  pulses  in  fiber  was  then  simulated  using  the  NLSE.  The  pulses  were 
then  recombined.  The  total  energy  transmitted  was  compared  to  the  energy  of  the  input 
pulse.  This  was  considered  to  be  the  transmission.  An  example  of  a  transmission  curve 
for  this  type  of  pulse  is  given  in  Figure  5.2.  This  curve  is  identical  with  the  results  of 
Duling  in  reference  78.  The  figure  displays  the  expected  dependence  of  the  energy 
needed  to  achieve  transmission  with  loop  length.  A  less  intense  pulse  needs  a  longer  loop 
to  build  up  a  phase  difference. 
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- 2  soliton  period  -  -  -  4  soliton  period - 6  soliton  period 

Figure  5.2  Pulse  transmission  as  a  function  of  pulse  energy  for  an  NOLM  with 
60/40  beam  splitter  and  no  birefringence.  The  squares  represent  the  predicted 
energies  required  for  peak  transmission  given  by  Equation  5.13.  The  figure  shows 

the  effect  of  changing  the  loop  length. 


The  picture  of  the  pulses  in  the  NOLM  forming  solitons  can  be  examined 
quantitatively  in  that  it  can  be  used  to  predict  the  location  of  transmission  maxima  in 
Figure  5.2  which  is  a  plot  of  the  transmission  curves  for  solitons  with  three  different 
lengths  of  the  NOLM.  This  can  be  done  by  using  the  following  simplified  picture  of  how 
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pulses  behave  in  the  NOLM.  The  pulses  will  be  transmitted  when  they  are  out  of  phase 
by  n  radians  when  they  return  to  the  beam  splitter.  I  derived  this  model  in  order  to  help 
understand  the  behavior  of  the  NOLM  and  also  to  enable  me  to  understand  how  changing 
the  parameters  of  the  NOLM  design  would  effect  their  behavior.  I  needed  some  way  of 
estimating  the  phase  shift  of  the  pulses  as  they  propagate  in  the  loop.  A  simple  way  of 
doing  this  is  to  assume  that  the  pulses  immediately  turn  into  solitons  with  the  same 
energy  as  the  counter-propagating  pulses.  I  can  then  assume  that  the  pulses  will 
experience  a  phase  shift  equal  to  that  the  solitons  would  experience.  The  assumption  that 
the  pulses  will  become  solitons  is  justifiable  since  it  is  observed  that  hyperbolic  secant 
shaped  pulses  with  amplitudes  0.5<A<1.5  will  evolve  into  solitons.4  For  the  model  to 
hold  both  of  the  counter-propagating  pulses  should  have  peak  intensities  in  this  range. 
This  will  be  true  if  the  amplitude  of  the  input  pulse  is  greater  than  (l-a)/2  and  less  than 
3 a/2.  For  a  60/40  beamsplitter,  for  instance,  the  range  of  input  intensities  will  be 
between  1.25  and  2.5.  The  assumption  that  these  pulses  immediately  behave  like  solitons 
is  an  oversimplification  and  this  explains  the  inaccuracies  of  the  model. 

We  must  now  relate  the  length  of  the  loop  to  a  phase  shift.  The  lengths  are  chosen 
in  units  of  the  soliton  period  z0.  The  phase  shift  (j)  of  a  soliton  is  given  by79 


where  A  is  the  pulse  intensity  amplitude  and  the  variable  \  is  distance  in  dispersion 
lengths.  As  mentioned  in  Equation  2.36  this  is  proportional  to  the  pulse  width.  It  is 
assumed  that  the  input  pulses  have  a  hyperbolic  secant  shape  and  an  energy  Ejn.  The 
width  of  the  input  pulses  are  held  constant  and  only  the  amplitude  is  varied.  These  are 
rapidly  reshaped  into  solitons  of  energy  Ep  for  which  the  amplitude  and  width  are 
coupled.  The  quantity  Ejn  is  given  by 
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Ein  =  2  J  pin  sec h2 [t]dt  =  2Pin .  (5.7) 

0 


In  contrast  the  energy  of  the  soliton  is  given  by 


Ep  =2  JPpSech' 


dt  =  2P  T  . 

P  P  : 


(5.8) 


where  Pp  and  Tp  are  the  pulse  height  and  width  of  the  soliton  formed.  In  the 
dimensionless  soliton  units  used  here  these  are  related  by 


If  we  use  Equations  5. 7-5. 9  and  assume  that  Ej„  =  Ep ,  we  get 

~  (5.10) 

1P 

and 

P„=(P„f.  (5-11) 

where  Tp  is  the  width  of  the  soliton  formed  by  the  input  pulse.  If  Equations  5.6,  5.10,  and 
5.11  are  used,  one  finds  that  the  phase  difference  between  the  two  counter-propagating 
pulses  becomes 

(2a  -  l)(P  ) 2 

(5.12) 

Solving  for  pulse  energy, 
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(5.13) 


P.„ 


2A(()  2 

(2a-%_ 


The  results  of  this  analysis  are  shown  as  square  markers  in  Figure  5.2.  It  shows  good 
agreement  between  the  model  and  the  results  of  the  simulations,  however,  the  model  does 
tend  to  underestimate  the  energy  required  for  switching.  This  is  to  be  expected  as  the 
pulse  will  have  had  less  time  to  be  reshaped  into  a  soliton. 

Figure  5.2  also  shows  that  transmission  is  highest  if  the  loop  has  a  length  of 
around  4  soliton  periods.  For  the  splitting  ratio  being  examined  in  this  figure,  a  loop  of 
this  length  will  have  maximum  transmission  for  input  pulses  with  an  energy  slightly 
higher  than  twice  that  of  a  unit  soliton.  This  means  that  the  pulses  formed  in  the  loop  by 
the  directional  coupler  are  already  close  to  being  unit  solitons  and  there  will  be  relatively 
little  reshaping  or  dispersive  wave  shedding  in  the  fiber  of  the  loop.  This  explains  the 
high  transmission. 

We  expect  that  the  effect  of  decreasing  the  splitting  ratio  will  be  to  require  a 
longer  loop  for  the  same  degree  of  switching.  This  is  due  to  the  relatively  small 
difference  in  the  energies  of  the  pulses  formed.  The  effect  of  making  this  change  is 
shown  in  the  Figure  5.3.  The  splitting  ratio  used  was  55/45.  The  lengths  chosen  were 
those  that  would  give  nearly  the  same  switching  lengths  as  in  Figure  5.2.  This  was 
determined  by  considering  Equation  5.13.  In  the  last  case  the  value  of  the  term  2a- 1  was 
0.2.  In  the  previous  case  it  was  0.1,  this  implies  that  one  may  keep  the  switching  lengths 
constant  by  doubling  the  loop  lengths.  The  predicted  locations  of  first  transmission  peaks 
by  the  soliton  model  are  shown  in  Figure  5.3  as  squares.  The  model  again  tends  to 
underestimate  the  energy  needed  for  switching. 
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Figure  5.3  Pulse  transmission  as  a  function  of  pulse  energy  for  an  NOLM  with  a 
55/45  beam  splitter  and  no  birefringence.  The  squares  represent  the  predicted 
energies  required  for  peak  transmission  given  by  Equation  5.13.  The  figure  shows 

the  effect  of  changing  the  loop  length. 
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Simulation  Results:  Birefringent  NOLM 


Previous  work  on  NOLMs  have  either  ignored  the  effect  of  birefringence  in  the 
fiber  or  have  examined  their  behavior  for  a  nearly  constant  signal.80,81’82  I  have  extended 
this  work  by  studying  pulse  propagation  in  a  birefringent  fiber  NOLM.  I  have  made  the 
idealization  of  fiber  with  a  uniform  birefringence.  In  particular,  since  fiber  lasers  have 
been  noted  to  be  very  sensitive  to  the  setting  of  the  polarization  controllers,  I  have  studied 
the  impact  of  the  settings  of  the  controllers.  This  effect  has  been  used  to  greatly  reduce 
the  pump  power  needed  to  start  a  figure  eight  laser.83  The  propagation  of  a  pulse  in  the 
anomalous  dispersion  regime  can  be  modeled  by  the  coupled  Nonlinear  Schrodinger 
Equations  (NLSEs). 

I  modeled  the  polarization  controller  in  a  fashion  described  by  Mortimore81.  The 
field  components  of  the  pulse  are  rotated  into  the  coordinate  frame  of  the  controller  with 
the  controller  axes  being  its  fast  and  slow  birefringent  axes.  A  phase  delay  is  imposed  on 
the  slow  axis  and  the  fields  are  rotated  back  into  the  original  coordinate  system.  The 
resulting  Jones  matrix  for  the  controller  is 


( E'  ^ 
vE'yy 


(3 

xx 
V  J  yx 


xy 


yy/ 


E, 

vEyy 


(5.14) 


Ex  and  Ey  are  the  electric  field  components  in  the  input  coordinate  frame.  The  coordinate 
system  used  here  was  that  the  slow  axis  of  the  fiber  was  the  x  coordinate  axis  and  the  fast 
axis  was  the  y  axis.  The  elements  of  the  transfer  matrix  are  given  by 

Jxx  =  e^sin2  0  +  cos2  9, 

Jxy  =  Jyx  =  (e‘*  -  l)  cos 0  sine,  (5.15) 

Jyy  =e'*  cos2  0  + sin2  0. 

The  variable  0  is  the  angle  between  the  slow  axis  of  the  fiber  and  the  slow  axis  of  the 
polarization  controller.  The  relative  phase  shift  caused  by  the  controller  on  its  slow  axis 
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when  compared  to  its  fast  axis  is  given  by  <j).  An  additional  subtlety  of  extreme 
importance  was  pointed  out  by  Mortimore.81  The  effect  of  the  loop  geometry  is  that  the 
components  of  the  electric  field  of  the  counter  propagating  pulses  in  the  plane  of  the  loop 
have  their  orientations  reversed  relative  to  each  other.  This  means  that  their  orientations 
are  not  the  same  when  they  enter  the  polarization  controller.  The  significance  of  this  will 
become  clear  when  the  results  of  the  simulations  are  examined. 

In  general,  all  angles  have  been  measured  relative  to  the  slow  birefringent  axis  of 
the  fiber.  There  are  three  other  angles  of  importance.  One  must  define  the  orientation  of 
the  input  pulse  and  the  angle  of  the  plane  of  the  loop.  As  mentioned  in  the  last  paragraph, 
the  component  of  the  light  in  the  plane  of  the  loop  will  have  its  orientation  reversed  in 
one  counter  propagating  pulse  relative  to  the  other.  Finally,  one  must  define  the  direction 
of  the  slow  axis  of  the  polarization  controller.  These  directions  are  shown  in  Figure  5.4. 

I  have  sought  to  find  an  NOLM  configuration  which  would  have  a  set  of  desirable 
properties.  Clearly,  the  device  needs  to  have  high  transmission,  but  the  loop  also  needs  to 
transmit  pulses  which  have  been  undistorted.  If  the  distortions  are  too  great  they  will 
cause  the  pulse  to  break  up  in  subsequent  transmission  through  fiber  or  in  a  subsequent 
transmission  through  another  NOLM.  In  addition  I  would  like  to  be  able  to  change  the 
transmission  of  the  NOLM  by  adjusting  the  polarization  controller.  This  is  an  attractive 
option  since  this  is  often  the  only  adjustment  to  the  NOLM  that  can  be  made  once  it  has 
been  assembled.  It  should  be  emphasized  that  almost  any  desired  Jones  matrix  can  be 
produced  within  a  polarization  controller.  For  this  reason,  one  can  assume  that  in  an 
experiment  one  could  adjust  the  input  polarization  into  some  fixed  state. 
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In  order  to  make  the  polarization  controller  function  in  this  way,  I  need  to  create 
an  NOLM  configuration  in  which  the  two  counter-propagating  pulses  are  treated  quite 
differently  by  the  controller.  This  situation  can  be  achieved  in  the  following  way.  One  of 
the  pulses  must  be  aligned  with  one  of  the  birefringent  axes  of  the  polarization  controller. 
There  must  be  an  angle  of  45  degrees  between  the  orientation  of  the  input  pulse  and  the 
plane  of  the  loop.  This  will  cause  a  90  degree  change  in  the  orientation  of  one  of  the 


orientation  of 

counterpropagating 

pulse 


input  orientation 


slow  axis 


plane  of  the  loop 


Figure  5.4  Relative  directions  of  important  planes 
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pulses  as  it  traverses  the  loop.  Hence  one  of  the  pulses  will  be  aligned  with  the  slow  axis 
of  the  controller,  and  the  other  with  the  fast  axis.  In  this  configuration,  a  simple  relative 
phase  shift  of  one  pulse  relative  to  the  other  can  be  produced. 

I  have  examined  situations  in  which  this  simple  phase  shift  is  produced.  Since  the 
polarization  controllers  can  be  changed  to  form  any  arbitrary  Jones  matrix,  the 
polarization  of  the  input  pulse  can  be  altered  to  any  arbitrary  state.  Figure  5.5  shows  an 
example  of  a  set  of  transmission  curves  where  the  pulse  energy  is  equally  divided 
between  the  two  birefringent  axes  of  the  fiber.  The  plane  of  the  loop  was  aligned  with  the 
slow  axis  of  the  fiber  as  was  the  slow  axis  of  the  polarization  controller.  The  group 
velocity  difference,  8,  was  set  equal  to  0.1.  The  pulse  energy  is  normalized  so  that  the 
energy  of  a  unit  soliton  is  unity  and  the  directional  coupler  had  a  60/40  splitting  ratio. 

The  input  pulse  had  a  form  (A/2)sech2[(l+2/3)1/2x],  where  A  is  the  amplitude  and  x  is 
the  normalized  time  for  a  unit  soliton.  This  is  a  form  clearly  related  to  that  in  Equation 
3.18.  The  decision  to  use  this  form  was  guided  by  the  need  for  the  pulses  to  be  at  least 
potentially  stable.  The  loop  was  two  soliton  periods  long,  for  the  same  unit  amplitude 
soliton.  The  figure  indicates  that  this  NOLM  configuration  can  control  the  loop’s 
switching  functions  in  the  fashion  which  I  desire.  When  the  curve  describing  pulse 
transmission  without  phase  delay  is  compared  to  that  when  there  is  a  90  degree  phase 
delay,  we  find  that  the  energy  needed  to  reach  maximum  transmission  has  been  cut  in  half 
and  low  energy  transmission  has  been  increased  by  more  than  ten-fold. 
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Figure  5.5  Pulse  transmission  as  a  function  of  pulse  energy  for  an  NOM  with  a 
60/40  beam  splitter.  The  fiber  has  a  birefringence  of  5=0.1,  and  the  loop  is  2  soliton 
periods  long.  The  squares  represent  the  predicted  energies  required  for  peak 
transmission  given  by  Equation  5.20.  The  figure  shows  the  effect  of  change  in  phase 

delay. 
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The  location  of  the  transmission  maximum  for  the  curve  will  be  examined  first. 
At  the  maximum  the  two  pulses  must  be  out  of  phase  by  %.  The  phase  shift  a  soliton 
experiences  as  it  propagates  is 


<l>  = 


%■ 


(5.16) 


The  energy  of  the  input  pulse  is  given  by 


E|n  =  2  J  Pin  sech2 
0 


(5.17) 


where  Pjn  is  the  pulse  height  of  the  input  pulse.  In  contrast  the  energy  of  the  soliton  is 
given  by 


(5.18) 


This  indicates  that  Equations  5.9-5. 1 1  hold  in  this  case  as  well.  If  these  are  combined 
with  Equation  5.18  one  finds  that  the  phase  difference  between  the  two  counter- 
propagating  pulses  becomes 


(5.19) 


where  £,i  is  the  length  of  the  loop  in  dispersion  lengths.  Solving  for  pulse  energy  yields 


(5.20) 
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The  results  of  this  simple  model  and  those  predicted  by  numerical  simulation  are  given  in 
Figure  5.5.  The  results  indicate  moderately  good  agreement,  especially  given  the  crude 
approximations  made. 

The  output  pulses  were  compared  to  the  form  of  the  input  pulses 


(5.21) 


where  To  is  the  altered  pulsewidth  in  normalized  units.  The  output  pulses  maintained 
nearly  this  shape  after  transmission  through  the  loop.  The  output  values  of  To  are  given  in 
Table  5.1.  We  see  that  the  effect  of  changing  the  phase  delay  is  to  make  the  pulses 
broader  as  well.  Since  the  pulses  have  less  energy  for  maximum  switching  when  a  phase 
delay  is  added,  the  pulses  will  tend  to  broaden  in  the  loop.  This  would  explain  the  trend 
we  observe.  The  NOLM  may  reshape  a  pulse  transmitted  through  it.85  For  an  NOLM 
without  birefringence,  this  reshaping  is  minimized  for  loops  that  are  four  soliton  periods 
long.86 


The  shift  in  phase  delay  produces  a  distortion  of  the  pulse.  First,  as  can  be  seen  in 
Figures  5. 6-5. 8  it  may  be  observed  that  the  output  pulses  are  nearly  hyperbolic  secant 
shaped  pulses.  The  solid  and  dashed  curves  represent  the  intensities  along  the  fast  and 
slow  axes.  The  wings  are  built  up  a  bit  as  the  phase  delay  is  increased.  The  peak 
intensity  is  decreased  as  the  phase  delay  goes  up,  but  this  is  simply  due  to  that  the  fact 
that  the  input  pulses  are  less  energetic.  This  may  be  an  additional  cause  of  the  pulse 
broadening  shown  in  Table  5.1.  This  result  has  been  reported  before  for  non-birefringent 
NOLMs.  The  pulses  appear  to  propagate  as  coupled  solitons  since  the  polarizations 
have  not  separated  in  time  in  spite  of  birefringence. 
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Table  5.1  Pulse  Width  as  a  function  of  phase  delay  for  a  vector  soliton  in  a  2  soliton 
period  long  NOLM  with  a  60/40  beam  splitter.  The  birefringence  is  5=0.1. 


phase  delay 

pulse  width 

0 

0.8 

7l/4 

0.9 

7t/2 

1.2 

A  second  feature  which  can  be  explained  analytically  is  the  minimum 
transmission.  It  can  be  found  because  the  phase  delay  is  only  caused  by  the  polarization 
controller.  Equation  5.4  tells  us  that 

Tmin  =  1  “  2a(!  -  «)[!  +  cosO)] .  (5.22) 

The  results  of  this  expression  and  those  found  from  computation  are  given  in  Table  5.2. 
The  results  show  excellent  agreement. 

Table  5.2  Numerical  and  analytical  low  intensity  transmission  for  a  2  soliton 
period  long  NOLM  with  a  60/40  beamsplitter.  The  birefringence  is  5=0.1.  The 
numerical  results  are  compared  to  the  predictions  of  equation  5.22. 


phase  delay 

numerical 

analytical 

disagreement 

0 

0.04 

0.04 

0% 

7l/4 

0.18 

0.18 

0% 

71/2 

0.50 

0.52 

4% 

69 


intensity 


-10 


-5 


0 

time 


5 


10 


slow  axis - fast  axis 


Figure  5.6  Output  pulse  shape  for  a  2  soliton  period  long  loop  mirror  without  phase 
delay.  The  birefringence  of  the  fiber  was  8=0.1. 


Figure  5.7  Output  pulse  shape  for  a  2  soliton  period  long  loop  mirror  rc/4  delay.  The 

birefringence  of  the  fiber  was  8=0.1. 
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Figure  5.8  Output  pulse  shape  for  a  2  soliton  period  long  loop  mirror  rc/2  delay.  The 

birefringence  of  the  fiber  was  §=0.1. 

A  second  consequence  of  pulse  reshaping  is  a  reduction  in  transmission  of  the 
pulse,  when  it  is  passed  through  an  NOLM  a  second  time.  This  has  been  examined 
experimentally  by  Nayar,  Finlayson,  and  Doran  .  They  demonstrated  that  the 
transmission  of  their  NOLM  was  about  70%  for  the  first  transmission  but  that  two 
concatenated  NOLMs  had  a  transmission  of  about  50%.  They  also  showed  that  the  pulses 
appeared  to  be  much  more  square  than  the  input  pulses.  Figure  5.9  shows  what  happens 
when  this  is  done  with  an  NOLM  with  the  configuration  that  was  discussed  above.  This 
phenomenon  is  a  useful  thing  to  study  since  multiple  transmissions  through  a  switch 
occur  in  a  fiber  communications  network,  or  in  a  laser.  In  both  cases,  the  pulses  would  be 
reshaped  by  various  devices  between  passes  through  the  NOLM,  so  this  is  merely  a 
starting  point  to  the  study  of  these  problems.  We  see  that  when  there  is  no  phase  delay, 
the  effect  of  the  second  pass  is  to  narrow  the  transmission  peak  and  to  reduce 
transmission  between  the  peaks.  For  (j)=7t/2  ,  however  the  transmission  peak  is 
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Figure  5.9  Two  pass  transmission  for  a  2  soliton  period  loop.  The  birefringence  of 

the  fiber  was  5=0.1 
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eliminated,  and  the  low  intensity  transmission  is  increased.  The  curve  for  <t>=7t/4  gives  an 
intermediate  result.  The  low  intensity  transmission  is  simply  the  square  of  the  one  pass 
low  intensity  transmission. 

The  explanation  for  the  behavior  of  the  first  transmission  peak  is  more  complex. 
This  is  due  to  the  fact  that  the  system  is  now  nonlinear.  The  transmission  peak  for  the 
pulse  which  does  not  experience  phase  delay  has  been  shifted  to  a  lower  energy  than  in 
the  one  pass  case.  It  is  located  at  5.2  rather  than  5.3  energy  units.  The  peak  transmission 
in  this  situation  is  0.978.  If  one  simply  assumes  that  the  reduction  is  linear,  one 
multiplies  the  starting  energy  by  the  transmission  for  that  energy  and  then  multiplies  the 
result  by  the  one  pass  transmission  of  this  intermediate  energy.  This  analysis  predicts  a 
two  pass  transmission  of  0.967  or  an  error  of  about  1%.  So  in  this  case  the  linear  analysis 
works.  If  a  7t/4  phase  delay  is  imposed,  the  analysis  breaks  down.  It  predicts  a 
transmission  of  0.893  which  is  far  higher  than  the  value  of  .597  predicted  by  the 
numerical  results.  This  lower  value  could  be  due  to  pulse  reshaping  The  pulse  width,  as 
seen  in  Table  5.1,  is  broader  than  in  the  case  without  phase  delay.  This  will  certainly  alter 
the  transmission  somewhat.  When  the  phase  delay  is  increased  to  tt/4  radians, 
transmission  is  greatly  reduced.  In  this  case,  the  pulse  is  still  broader  and  can  be  seen  to 
be  somewhat  distorted. 
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6  soliton  period 


Figure  5.10  Pulse  transmission  as  a  function  of  phase  delay  in  an  NOLM  without  a 
polarization  controller.  The  birefringence  of  the  fiber  was  8=0.1.  The  NOLM  had  a 
60/40  beam  splitter,  and  a  vector  soliton  input  pulse.  The  figure  shows  the  effect  of 
changing  the  loop  length.  The  squares  show  the  pulse  energies  for  peak 
transmission  given  by  equation  5.20. 


The  expected  result  of  making  the  loop  longer  is  to  reduce  the  energy  required  for 
transmission.  Figure  5.10  shows  the  result  of  carrying  out  the  simulation  of  the  previous 
case,  by  varying  the  loop  length  using  no  phase  delay  due  to  the  controller.  The  low 
intensity  transmission  is  almost  exactly  what  is  predicted  analytically,  as  shown  in  Table 
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5.2  The  results  again  show  good  agreement.  The  comparison  of  the  expected  peak 
location  to  the  results  of  the  simple  model  improve  as  the  loop  becomes  longer.  This  may 
be  attributed  to  the  fact  that  the  pulse  is  reshaped  into  a  soliton  soon  after  entering  the 
loop  and  propagates  in  this  form  for  a  longer  portion  of  the  total  distance.  Under  these 
conditions,  the  assumptions  of  the  simple  model  are  more  nearly  correct. 

The  width  of  the  output  pulse  is  dependent  on  the  length  of  the  loop.  Table  5.3 
shows  this  effect.  The  output  pulse  becomes  broader  as  the  length  of  the  loop  increases. 
This  is  similar  to  the  results  in  the  non-birefringent  case.  The  pulses  are  less  energetic  for 
the  longer  loops,  hence  they  spread  in  the  loops.  This  results  in  longer  output  pulses. 

Table  5.3  Pulse  width  as  a  function  of  loop  length  for  a  vector  soliton.  The  loop  has 
a  60/40  beam  splitter.  The  birefringence  is  8=0.5. 


loop  length 

To 

2 

0.8 

4 

1.1 

6 

1.3 
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Figure  5.11  Pulse  transmission  as  a  function  of  pulse  energy  for  a  2  soliton  period 
long  loop.  The  NOLM  had  a  birefringence  of  8=0.5  and  a  vector  soliton  input.  The 
figure  shows  the  effect  of  changing  the  phase  delay.  The  squares  show  the  pulse 
energies  for  peak  transmission  given  by  equation  5.20. 
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Table  5.4  Low  intensity  transmission  for  a  2  soliton  period  long  NOLM 
with  a  60/40  beamsplitter.  The  birefringence  was  5=0.5.  The  numerical 
results  are  compared  to  the  predictions  of  equation  5.22. 


phase  delay 

numerical 

analytical 

disagreement 

0 

0.04 

0.04 

0% 

tc/4 

0.12 

0.18 

50% 

tc/2 

0.32 

0.52 

63% 

The  effect  of  increasing  birefringence,  which  is  described  by  the  parameter  5,  has 
been  examined.  In  general,  the  effect  of  increased  birefringence  is  to  reduce  the 
functionality  of  the  NOLM  switch.  This  is  due  to  the  tendency  of  birefringence  to  distort 
the  vector  solitons,  making  them  less  stable.  The  primary  form  of  distortion  is  that  the 
light  in  the  two  polarizations  separate;that  is  to  say,  the  sub-pulses  along  the  two 
birefringent  axes  tend  to  separate  from  each  other.  This  tendency  is  reduced  by  the 
soliton  trapping  effect  which  causes  light  in  the  two  polarizations  to  couple.  The  degree 
of  separation  is  dependent  upon  the  birefringence  and  on  the  energy  in  the  soliton.  This 
effect  was  described  in  the  chapter  on  birefringence.  This  effect  has  a  variety  of 
manifestations,  as  will  be  seen.  Figure  5.11  shows  the  result  of  making  5=0.5  and 
examining  the  effect  of  changing  the  phase  delay.  The  first  effect  to  be  noted  is  that  the 
low  intensity  transmission  no  longer  agrees  with  the  cw  analysis.  This  is  shown  in  Table 
5.4,  which  shows  that  the  pulse  width  is  relatively  insensitive  to  changes  in  phase  delay 
for  this  case.  In  addition, we  can  see  that  the  predictions  of  the  location  of  the 
transmission  peaks  become  far  worse.  These  predictions  are  indicated  in  Figure  5.11. 
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Table  5.5  Pulse  width  as  a  function  of  phase  delay  for  a  2 
soliton  period  long  NOLM  with  a  60/40  beamsplitter.  The 
briefringence  is  8=0.5. 


phase  delay 

pulse  width 

0 

1.1 

7l/4 

1.0 

7t/2 

1.0 

' 

The  explanation  for  the  reduction  in  low  intensity  transmission  can  be  found  by 
examining  how  the  polarization  controller  affects  pulse  shape.  Figures  5.12-5.14  show 
the  case  of  a  hyperbolic  secant  pulse  in  the  low  intensity  case.  The  dashed  and  solid  lines 
represent  intensity  along  the  two  birefringent  axes.  The  loop  is  two  soliton  periods  long. 
The  polarization  controller  is  on  one  end  of  the  loop.  Of  the  two  pulse  traveling  in  the 
loop  mirror,  one  will  encounter  the  controller  first,  and  then  be  modified  by  the  fiber  in 
the  loop.  The  other  will  be  altered  by  the  polarization  controller  only  after  passing 
through  the  loop.  The  pulse  which  first  is  passed  through  the  polarization  controller  and 
then  passes  through  the  loop  is  undistorted.  The  pulse  which  passes  through  the  loop 
first,  however  is  dramatically  distorted  by  the  polarization  controller.  Figure  5.12  shows 
the  pulse  after  it  has  passed  through  the  loop  but  before  it  is  affected  by  the  polarization 
controller.  It  has  been  broadened  by  linear  dispersion  and  the  two  polarizations  are 
separating,  but  it  still  has  a  hyperbolic  secant  shape.  This  is  the  same  as  the  counter 
propagating  pulse  which  first  passes  through  the  polarization  controller  and  then  through 
the  fiber.  Figure  5.13  shows  the  the  pulse  which  has  been  passed  through  the  fiber 
through  the  controller.  Large  satellite  pulses  have  been  produced.  The  origin  of  these 
satellites  can  be  found  in  Equations  5.14  and  5.15,  which  describe  the  polarization 
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controller.  In  the  case  where  the  pulse  energy  is  equal  along  both  axes  (9=7r/4),  the 
polarization  controller  will  not  change  the  polarization  orientation  of  the  pulse  but  simply 
acts  to  impose  a  phase  delay  on  one  axis.  Figure  5.12  shows  however,  that  the  pulse  does 
not  in  general  have  equal  intensities  along  both  axes  at  any  given  moment.  Under  this 
condition,  the  controller  couples  light  between  the  two  input  polarizations.  This  explains 
the  distorted  pulse.  The  distortion  in  the  output  pulse  is  shown  in  Figure  5.14. 


Figure  5.12  Shape  of  a  low  energy  pulse  in  the  loop.  This  pulse  was  affected  by  the 
polarization  controller  at  the  beginning  of  its  passage  through  the  fiber  loop.  The 
birefringence  of  the  fiber  was  5=0.5.  There  is  a  phase  delay  of  tc/4.  The  loop  is  two 

soliton  periods  long. 
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Figure  5.13  Shape  of  a  low  energy  pulse  in  the  loop.  This  pulse  was  affected  by  the 
polarization  controller  at  the  end  of  its  passage  through  the  fiber  loop.  The 
birefringence  of  the  fiber  was  8=0.5.  There  is  a  phase  delay  of  u/4.  The  loop  is  two 

soliton  periods  long. 
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Figure  5.14  Shape  of  a  low  energy  pulse  transmitted  by  an  NOLM.  The 
birefringence  of  the  fiber  was  5=0.5.  There  is  a  phase  delay  of  7r/4.  The  loop  is  two 

soliton  periods  long. 


The  distortion  of  one  of  the  counter  propagating  pulses  causes  the  change  in 
transmission.  The  relationship  between  the  two  can  be  seen  by  taking  the  magnitude  of 
both  sides  of  Equation  5.3  to  obtain 


|Eoutf  =a|E3f  +(l-a)|E4|2-2Re  ia^  (l-a) 


E3E4e11 


(5.21) 


This  equation  expresses  the  output  of  the  NOLM,  where  the  last  term  describes  the  effect 
of  interference  of  the  two  counterpropagating  pulses.  If  the  magnitude  of  the  pulses  can 
be  related  to  the  input  fields  by  Equation  5.1,  then  Equation  5.21  reduces  to  Equation  5.4. 
But  since  one  of  the  pulses  is  distorted  by  the  polarization  controller,  the  instantaneous 
transmission  is  altered.  The  low  intensity  transmission  cited  in  Table  5.4,  represents  the 
averaging  of  the  instantaneous  transmission  over  time.  Indeed  if  one  uses  Equation  5.1  in 
the  last  equation,  but  assumes  that  coherence  between  the  two  pulses  has  been  eliminated, 
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one  obtains  a  transmission  of  0.52  if  a=0.6.  This  is  what  we  see  in  Table  5.4  for  the  case 
where  <j>=7t/2.  So  in  this  case  the  interference  term  is  making  no  net  contribution  to  the 
pulse  transmission. 

The  reason  that  this  phenomenon  did  not  appear  for  the  case  of  8=0. 1 ,  can  be  seen 
in  Figures  5.15-5.17.  The  dashed  and  solid  lines  represent  intensity  along  the  two 
birefringent  axes  and  the  loop  is  two  soliton  periods  long.  The  NOLM  has  a  60/40 
beamsplitter.  These  show  the  pulses  under  the  same  conditions  except  with  this  lower 
value  of  8.  We  can  clearly  see  here  that  under  these  conditions  the  light  in  the  two 
polarizations  separate  less.  The  pulses  in  Figure  5.15  display  only  about  one  fifth  the 
separation  as  in  Figure  5.12.  This  is  to  be  expected  from  the  difference  in  birefringence. 
Soliton  trapping  does  not  occur  due  to  the  low  pulse  intensities.  The  much  smaller  pulse 
separation  means  that  Equation  5.22  is  a  much  better  approximation  and  the  simple 
theory  holds  more  exactly.  The  effect  of  the  controller  can  be  seen  in  Figure  5.16.  It  is 
manifested  in  an  increase  in  pulse  separation  rather  than  by  the  large  wings  shown  in 
Figure  5.13,  which  is  the  higher  birefringence  case.  The  final  output  pulse  is  shown  in 
Figure  5.17.  The  pulse  appears  undistorted  in  the  8=0.1  case.  The  8=0.5  results,  shown 
in  Figure  5.14,  again  show  large  wings. 
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Figure  5.15  Shape  of  a  low  energy  pulse  in  the  loop.  This  pulse  was  affected  by  the 
polarization  controller  at  the  beginning  of  its  passage  through  the  fiber  loop.  The 
birefringence  of  the  fiber  was  8=0.1.  There  is  a  phase  delay  of  rc/4.  The  loop  is  two 

soliton  periods  long. 
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Figure  5.16  Shape  of  a  low  energy  pulse  in  the  loop.  This  pulse  was  affected  by  the 
polarization  controller  at  the  end  of  its  passage  through  the  fiber  loop.  The 
birefringence  of  the  fiber  was  5=0. 1 .  There  is  a  phase  delay  of  tc/4.  The  loop  is  two 

soliton  periods  long. 
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Figure  5.17  Shape  of  a  low  energy  pulse  transmitted  by  an  NOLM.  The 
birefringence  of  the  fiber  was  8=0. 1 .  There  is  a  phase  delay  of  td/4.  The  loop  is  two 

soliton  periods  long. 


The  same  pulse  reshaping  phenomenon  helps  to  explain  the  fact  that  the  peak 
transmission  for  the  8=0.5  case  declines  as  the  phase  delay  increases.  It  also  manifests 
itself  by  the  altering  of  the  NOLM  pulse  output  shape.  This  is  seen  in  Figures  5.18-5.20. 
These  figures  show  the  output  pulses  transmitted  by  the  NOLM  as  the  phase  delay  is 
increased.  The  loop  is  two  soliton  periods  long  and  it  has  a  60/40  beam  splitter.  The  solid 
and  dashed  lines  indicate  intensity  along  the  two  birefringent  axes.  These  show  the  effect 
of  adding  more  phase  delay.  The  difference  between  the  results  shown  here  and  those 
seen  for  low  intensity  pulses  is  that  the  pulses  in  the  two  polarizations  will  couple.  This 
nonlinear  effect  will  reduce  the  separation  and  hence  the  distortion  of  the  pulse  by  the 
polarization  controller.  At  first  there  is  simply  some  separation  of  the  polarizations. 
When  the  delay  is  7t/4,  the  pulse  has  a  hyperbolic  secant  shape  which  has  reduced 
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polarization  separation.  This  shows  clearly  how  the  soliton  dragging  effect  helps  the 
NOLM  function  in  this  case.  In  the  7i/2  case,  however,  there  is  increased  separation  and 
satellite  pulses  develop.  The  satellites  are  still  much  smaller,  relatively,  than  in  the  low 
intensity  case. 


Figure  5.18  Output  pulse  shape  for  a  2  soliton  period  long  loop  mirror.  There  is  no 
phase  delay.  The  birefringence  of  the  fiber  was  8=0.5. 
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Figure  5.19  Output  pulse  shape  for  a  2  soliton  period  long  loop  mirror.  The 
birefringence  of  the  fiber  was  8=0.5.  There  is  a  phase  delay  of  7t/4. 
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Figure  5.20  Output  pulse  shape  for  a  2  soliton  period  long  loop  mirror.  The 
birefringence  of  the  fiber  was  8=0.5.  There  is  a  phase  delay  of  Till. 

The  fact  that  the  polarization  controller  distorts  the  pulse  can  also  be  clearly  seen 
in  the  way  it  reduces  the  multiple  pass  transmission,  as  seen  in  Figure  5.21 .  In 
comparison  to  Figure  5.9,  the  transmission  peak  for  a  tc/4  phase  delay  has  been  almost 
eliminated.  The  low  intensity  transmission  has  increased,  as  is  expected.  The  peak 
transmission  for  the  switch  without  phase  delay  has  been  reduced  due  to  pulse  distortion. 
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Figure  5.21  Two  pass  pulse  transmission  for  a  vector  soliton  in  a  2  soliton  period 
loop.  The  birefringence  of  the  fiber  was  8=0.5. 

Figure  5.22  shows  the  effect  of  changing  the  loop  length  for  the  case  of  an  NOLM 
with  8=0.5.  The  behavior  of  the  NOLM  is  much  the  same  as  in  the  other  cases.  The 
effect  of  the  birefringence  is  to  increase  the  energy  required  to  switch  pulses.  The  peak 
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transmissions  are  reduced  due  to  pulse  distortion.  This  also  appears  in  Table  5.6  which 
shows  how  the  pulse  widths  increase  as  the  loop  becomes  longer.  Comparison  with 
Table  5.3  shows  that  this  effect  is  greater  than  in  the  case  with  less  birefringence.  The 
increased  birefringence  causes  the  pulses  along  the  two  birefringent  axes  to  separate  more 
in  time.  For  this  reason,  when  they  are  recombined  by  the  polarization  controller  the 
output  pulse  will  be  wider.  The  increase  in  birefringence  also  means  that  the  pulses  will 
propagate  like  solitons  for  a  reduced  part  of  their  passage  through  the  loop.  The  degree  to 
which  the  pulses  in  the  loop  do  not  behave  like  true  solitons  contributes  to  the  breakdown 
in  the  accuracy  of  the  soliton  model  for  predicting  pulse  energies  for  switching,  as  seen  in 
Figure  5.22.  It  should  be  noted  that  the  energies  needed  for  peak  transmission  are 
increased  from  the  5=0.1  case.  This  indicates  that  the  phase  difference  builds  up  more 
slowly  when  the  birefringence  increases.  Since  the  birefringence  parameter  does  not 
appear  in  Equation  5.20,  which  predicts  the  transmission  energy,  we  must  attribute  the 
change  to  increased  non-soliton  like  propagation. 
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Figure  5.22  Transmission  of  a  vector  soliton  as  a  function  of  soliton  energy.  The 
birefringence  of  the  fiber  was  8=0.5.  The  figure  shows  the  effect  of  changing  the 
loop  length.  The  squares  represent  the  predicted  energies  required  for  peak 
transmission  given  by  Equation  5.20. 
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Table  5.6  Pulse  width  as  a  function  of  loop  length  soliton  for  an  NOLM  with 
a  60/40  beamsplitter.  The  birefringence  is  8=0.5. 


loop  length 

To 

2 

1.1 

4 

1.4 

6 

1.8 

The  configuration  that  has  been  examined,  with  the  input  pulse  energy  evenly 
divided  between  the  two  birefringent  axes,  is  not  the  only  one  that  displays  switching. 

The  general  requirement  is  that  the  input  pulse  be  at  a  7i/4  angle  to  the  plane  of  the  loop 
and  that  one  of  the  axes  of  the  polarization  controller  be  aligned  with  the  input  pulse. 
Other  configurations  will  be  examined  much  more  briefly.  Figure  5.23  shows  the  case  in 
which  the  input  pulse  is  aligned  with  the  slow  axis.  The  central  change  in  this 
configuration,  compared  to  the  last  is  that  the  counter-propagating  pulses  are  now  moving 
with  different  group  velocities.  For  this  reason,  the  interference  between  them  is  reduced 
when  they  recombine.  This  is  manifested  in  the  reduced  peak  transmissions  observed  in 
Figure  5.23  when  compared  to  Figure  5.5.  This  effect  also  appears  in  the  increased  low 
intensity  transmission  displayed.  Figure  5.24  shows  the  two-pass  transmission  when 
8=0. 1 .  It  indicates  that  this  configuration  is  less  sensitive  to  pulse  distortion  in  the 
polarization  controller  than  the  case  in  which  the  light  is  evenly  divided  between  the 
birefringent  axes,  especially  in  the  case  where  the  phase  delay  is  7i/2.  The  pulse  distortion 
in  the  vector  soliton  case  was  caused  by  the  light  of  different  polarization  separating.  In 
this  case,  however,  this  does  not  happen,  so  the  polarization  controller  can  not  distort  the 
pulse.  This  explains  why  the  one  pass  transmission  is  more  uniform.  The  energy 
required  for  peak  two-pass  transmission  is  not  greatly  altered  by  the  polarization 
controller,  so  it  may  not  be  a  useful  adjustment. 
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The  effect  of  increasing  birefringence  is  far  more  dramatic.  Figure  5.25  displays 
this  effect.  Here  the  increased  pulse  separation  greatly  reduces  pulse  transmission.  This 
effect  can  be  understood  in  terms  of  the  ideas  given  above.  The  pulses  separate,  and 
hence  they  do  not  interfere  as  much.  This  would  lead  to  the  dramatic  reduction  in  peak 
transmission  shown  in  the  figure.  Figure  5.25  also  indicates  that  the  two  pass 
transmission  is  insensitive  to  changing  the  phase  delay  as  is  the  low  intensity 
transmission.  Hence  the  utility  of  adjusting  the  polarization  controller  is  lost  in  this  case. 
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Figure  5.23  NOLM  pulse  transmission  for  an  input  pulse  and  polarization 
controller  aligned  to  the  slow  axis  of  the  fiber.  The  plane  of  the  loop  is  tilted  by  7r/4 
from  this  and  the  birefringence  of  the  fiber  was  5=0.1.  The  loop  is  four  soliton 
periods  long.  The  figure  shows  the  effect  of  change  the  phase  delay. 
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Figure  5.24  Two  pass  NOLM  pulse  transmission  for  an  input  pulse  and  polarization 
controller  aligned  to  the  slow  axis  of  the  fiber.  The  plane  of  the  loop  is  tilted  by  tc/4 
from  this  and  the  birefringence  of  the  fiber  was  5=0.1.  The  loop  is  four  soliton 

periods  long. 
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Figure  5.25  Two  pass  NOLM  pulse  transmission  for  an  input  pulse  and  polarization 
controller  aligned  to  the  slow  axis  of  the  fiber.  The  plane  of  the  loop  is  tilted  by  7t/4 
from  this  and  the  birefringence  of  the  fiber  was  6=0.5.  The  loop  is  four  soliton 

periods  long. 


Transmission  curves  for  an  arbitrary  input  polarization  can  be  expected  to  be  more 
complex  since  the  pulses  will  experience  polarization  rotation  as  well  as  separation  of 
light  in  different  polarizations.  Figure  5.25  shows  the  polarization  curves  for  the  situation 
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in  which  the  input  pulse  is  at  an  angle  of  tt/8  to  the  slow  axis  of  the  fiber.  The  slow  axis 
of  the  controller  is  aligned  with  the  pulse  and  the  plane  of  the  loop  is  at  an  angle  of  -71/8. 
Pulses  with  this  input  angle  will  be  subject  to  intensity  dependent  polarization  rotation. 
Their  behavior  can  be  expected  to  be  quite  complex.  In  addition,  their  pulses  will  be 
subject  to  distortion  by  the  polarization  controller.  We  observe  that  the  peak  transmission 
is  reduced.  This  may  be  ascribed  to  reduce  interference  at  the  beamsplitter  due  to  the 
alteration  of  the  polarization  in  the  loop  and  to  the  separation  of  the  light  along  the  two 
different  axes.  The  separation  of  light  along  the  two  axes  manifests  itself  in  increased 
low  intensity  transmission,  compared  with  the  cw  case. 
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Figure  5.26  NOLM  pulse  transmission  for  an  input  pulse  and  polarization 
controller  at  a  71/8  angle  from  the  slow  axis  of  the  fiber.  The  plane  of  the  loop  is 
tilted  by  -tc/8  from  this  and  the  birefringence  of  the  fiber  was  8=0. 1 .  The  loop  is  four 

soliton  periods  long. 


Figure  eight  lasers  are  mode-locked  using  nonlinear  optical  loop  mirrors.  These 
fiber  components  have  also  been  used  for  a  number  of  other  applications  as  discussed  in 
the  beginning  of  this  chapter.  The  operation  of  the  NOLM  was  simulated  and  compared 
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to  previously  published  computational  results  to  verify  the  models  validity.  A  novel 
model  was  developed  to  predict  the  input  pulse  energy,  and  the  maximum  amount  of  the 
input  energy  to  be  transmitted  by  the  loop.  I  then  extended  these  results  by  examining  the 
effects  of  birefringence  on  the  functioning  of  the  loop.  It  was  found  that  the  transmission 
of  low  energy  pulses  could  be  increased  by  changing  the  relative  phase  delay  between  the 
two  counter-propagating  pulses  in  the  loop.  This  effect  could  be  explained  quantitatively 
as  long  as  the  birefringence  was  not  too  large.  The  break  down  of  the  model  was 
explained.  In  addition,  a  model  was  developed  to  predict  peak  pulse  transmission 
energies,  based  on  vector  soliton  propagation.  The  distortion  of  pulses  by  the  NOLM 
and  the  effect  of  this  on  further  transmissions  through  the  loop  was  examined.  These 
effects  have  an  influence  on  the  behavior  of  the  figure  eight  laser.  Hence  a  good 
foundation  has  been  established  for  the  simulation  of  this  system. 
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Chapter  Six  :  Figure  Eight  Fiber  Lasers 

The  components  described  in  previous  sections  can  be  combined  to  form  a  fiber 
laser.  This  chapter  starts  with  a  short  review  of  previous  work  on  fiber  lasers.  The  theory 
of  sideband  formation  is  reviewed  in  greater  detail.  This  is  the  most  successful  theory  in 
this  field  and  it  has  been  used  to  help  design  other  lasers.  I  then  go  on  to  simulate  a  fiber 
laser  mode-locked  by  an  NOLM  switch.  I  show  that  I  can  reproduce  in  the  simulation 
many  of  the  phenomena  observed  experimentally.  With  insights  gained  by  this  task  I 
have  proposed  the  concept  of  a  dispersion  balance  figure  eight  laser  which  overcomes 
some  of  the  limitations  of  previous  figure  eight  designs.  I  have  extensively  simulated  a 
laser  of  this  type  and  the  concept  appears  to  work  well.  I  then  compare  the  results  of  this 
work  to  the  sideband  formation  theory  described  before. 


Previous  Work 

A  laser  which  uses  the  NOLM  switch  is  shown  in  Figure  5.1.  Such  a  laser  was 
first  built  by  Avramopoulos  et  al.88  Another  laser  of  the  same  type  was  recently  reported 
by  Wu.  A  similar  type  of  laser,  based  on  the  NALM  switch,  has  also  been  built,  and 
much  more  work  has  been  done  on  these.  The  NALM  lasers  were  first  demonstrated  by 
Duling12,90,  and  by  Richardson  et  al.91  Pulses  as  short  as  100  fs  have  been  produced  in 
these  systems.92 

These  lasers  can  have  very  erratic  pulse  repetition  rates.93  This  is  one  of  the  most 
significant  problems  in  making  them  practical.  Others  have  found  that  the  NALM  laser 
has  several  distinct  repetition  rate  regimes  which  change  as  pump  power  is  increased.  At 
first,  the  laser  will  support  one  pulse  in  the  loop  which  means  that  the  output  repetition 
rate  of  the  laser  will  be  the  round  trip  time.  If  the  pump  power  is  increased,  one  observes 
a  chaotic  pulse  output  rate.  This  is  characteristic  of  harmonically  mode-locked  lasers.94 
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If  the  pump  power  is  increased  even  more,  bursts  of  chaotically  spaced  pulses  are 
observed.  It  will  be  shown  that  I  can  explain  all  three  operating  regimes. 

A  number  of  means  have  been  used  to  stabilize  the  repetition  rate.  A  fiber  sub¬ 
ring  can  be  introduced  into  the  cavity.95,96  The  sub-ring  acts  like  a  Fabry-Perot  etalon  to 
stabilize  the  modes  of  the  laser.  This  is  similar  to  the  means  used  to  stabilize  actively 
mode-locked  fiber  lasers.97  It  has  also  been  found  that  introducing  a  small  amount  of 
extracavity  feedback  can  help  to  stabilize  a  figure  eight  laser.98  Fiber  lasers  in  a  ring 
geometry  have  been  stabilized  with  phase  modulation.99  A  ring  laser  was  stabilized  by 
adjustments  made  to  polarization  controllers  in  the  loop.100  The  author  attributed  this 
stabilization  to  a  long  range  repulsive  interaction  between  the  solitons. 


polarization 


Figure  6.1  NOLM  Figure  Eight  Laser 

In  addition  to  mode-locking  with  loop  mirrors,  passive  mode  locking  in  fiber 
lasers  has  been  demonstrated  by  means  polarization  rotation  filtering.  This  has  been 
shown  in  a  Fabry-Perot  geometry.13,101  Polarization  rotation  for  mode-locking  has  been 
extensively  studied  in  a  ring  geometry  and  has  been  used  to  develop  a  type  of  laser  known 
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as  the  stretched  pulse  ring  laser.14  Nonlinear  switching  has  been  demonstrated  in  dual¬ 
core  optical  fibers.10  This  has  been  proposed  as  a  means  of  mode-locking  a  laser  as 
well.15 

The  model  developed  here  involves  the  simulation  of  the  laser  by  modeling  the 
propagation  of  a  pulse  through  optical  fiber.  Little  work  has  been  done  to  simulate  figure 
eight  lasers  in  this  way.  Bulushev  et  al.  102  simulated  a  figure  eight  laser  mode-locked  by 
an  NOLM  loop.  They  modeled  pulse  propagation  in  the  loop  but  not  in  the  amplifier, 
which  they  treated  as  a  lumped  gain.  It  will  be  shown  that  pulse  shaping  in  the  amplifier 
is  quite  important.  They  also  assumed  an  unrealistically  high  value  of  saturation  energy. 
This  allowed  them  to  demonstrate  self-starting.  Specifically,  they  assumed  an  amplifier  .* 
of  the  form 


ANo(cd) 


where  their  variable  £  corresponds  to  the  saturation  energy.  They  used  a  value  of  1 .8.  As 
was  shown  in  chapter  4,  even  when  using  a  300  fs  pulse  as  the  basis  of  the  normalized 
units,  the  saturation  energy  should  be  about  10'6.  Hence  their  results  are  invalid. 

Some  numerical  results  were  reported  by  Stentz  and  Boyd.103  They  did  not 
describe  their  model  except  to  say  that  it  involved  birefringent  pulse  propagation.  They 
assumed  that  the  birefringent  axes  of  the  fiber  rotated  at  random  intervals  in  the  fiber. 
They  demonstrated  that  this  eliminated  the  effect  of  nonlinear  polarization  rotation  in  the 
fiber. 


The  NALM  laser  was  modeled  by  Tzelpis  et  al.104  They  included  the  effect  of 
propagation  in  the  amplifier,  however,  they  also  used  a  very  high  value  of  saturation 
energy.  This  renders  their  results  suspect  as  well. 
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A  model  was  described  by  Chen  et  al.105  which  used  pulse  propagation  through  a 
stretch  of  birefringent  fiber  to  show  the  effect  of  pulse  rotation  in  an  actively  mode- 
locked  laser.  They  use  a  lumped  amplifier  model  along  with  a  lumped  saturable  absorber. 
They  incorporated  a  Fabry-Perot  etalon  in  order  to  stabilize  the  pulse  repetition  rate. 

They  were  able  to  foresee  the  need  to  incorporate  active  mode-locking  into  this  laser  in 
spite  of  the  saturable  absorption  used. 

Theoretical  and  experimental  results  have  been  reported  by  Richardson  et  al.106 
They  studied  the  result  of  amplifying  pulses  after  they  have  been  emitted  by  the  laser. 

They  demonstrated  that  when  the  total  gain  of  the  amplifier  is  increased,  one  can  obtain 
short  pulses  due  to  pulse  compression  in  the  amplifier. 

An  alternative  method  is  to  study  the  laser  by  means  of  an  NLSE  with  cyclical 
perturbations. 107,108,1 09,1 10,1 11,1 12  These  papers  show  how  the  periodic  amplification  of  the 
pulse  can  generate  spectral  sidebands  and  study  the  role  of  dispersion  in  breaking  up 
pulses.  The  spectral  side  bands  are  built  up  by  constructive  interference  between  the 
soliton-like  pulse  in  the  laser  and  the  cw  background.  The  background  is  built  up  by  the 
process  of  dispersive  wave  shedding.  When  the  pulse  is  perturbed  from  being  a  soliton  if 
the  perturbation  is  small,  the  pulse  will  be  reshaped  back  into  a  soliton.  The  process  will 
result  in  some  light  being  shed  by  the  pulse  as  cw  background.  The  amount  of  light  shed 
will  be  proportional  to  the  size  of  the  perturbation. 

The  location  of  the  sidebands  can  be  located  analytically  for  ring  lasers.  This 
analysis  will  be  adapted  to  figure  eight  lasers,  so  it  is  included  here.  The  frequencies  at 
which  the  dispersive  wave  radiation  interferes  constructively  with  the  soliton  must  be 
found.  The  phase  shift  of  the  dispersive  wave  radiation  will  be  (3(co)L  where  L  is  the 
length  of  the  ring.  The  soliton  phase  shift  is  %\H  as  discussed  above.  The  variable  %\  is 
the  length  of  the  cavity  in  dispersion  lengths.  The  soliton  envelope  function  only  contains 
the  slowly  varying  information,  which  is  dominated  by  the  second  order  term  of  the 
expansion  of  the  dispersion.  I  must  add  other  components  which  have  been  left  out.  The 
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actual  soliton  field  phase  will  contain  the  factor  of  exp[-i(3(co0)L],  If  this  is  expanded  to 
first  order,  the  phase  difference  becomes 


A(j)(co  +  6co)  =  P0  +  Pj  (Soo)  + - P(co  +  8co)  L . 

2Ld 


where  8co  -  co-co0  and  Ld  is  the  dispersion  length.  For  constructive  interference  we  must 
have 

A(f >  =  2m7t ,  (6.2) 

where  m  is  some  integer.  If  the  term  P(cd+8cd)  is  expanded  as  a  Taylor  series  and  inserted 
in  the  last  equation,  we  obtain 


L  =  2m7t . 


(6.3) 


Solving  for  §0,  the  result  is 


The  soliton  pulse  width  is  Tq2  —  Ld  P2,  so  the  last  equation  can  be  expressed  as 
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(  4m7iLD  )  ^ 

l  iT^_1J 


=  8Q, 


(6.5) 


where  5Q  is  in  normalized  units. 

The  analysis  done  by  others  has  shown  that,  in  general,  stable  operation  with  a 
fiber  laser  can  not  be  achieved  if 

L  <  4tiLd  (6.6) 

This  is  equivalent  to  saying  that  the  spacing  of  the  sidebands  can  not  be  allowed  to  be 
reduced  to  zero,  from  Equation  6.5.  That  is,  the  phase  change  for  the  soliton  must  be  less 
than  27r.  I  have  found  that  this  seems  to  hold  for  situations  in  which  perturbations  are 
producing  a  moderate  amount  of  dispersive  wave  radiation,  however,  if  the  perturbations 
are  reduced,  the  limit  can  be  eliminated. 

A  model  known  as  additive  pulse  mode  locking  has  been  the  basis  of  a  great  deal 
of  research.  Figure  eight  lasers  are  an  example  of  additive  pulse  mode  locking  (APM) 
and  have  been  analyzed  in  this  fashion.  The  vast  majority  of  work  involving  additive 
pulse  mode-locking  has  addressed  ring  lasers.  APM  theory  involves  constructing  a 
master  equation  of  the  laser  and  finding  solutions  to  it.  The  master  equation  is 
constructed  by  assuming  that  the  pulse  is  subject  to  only  small  perturbations  as  it  goes 
through  the  loop.  This  is  not  the  case  for  figure  eight  lasers  since  the  pulse  in  the 
amplifier  will  have  approximately  twice  the  amplitude  of  the  pulse  in  the  loop  mirror.  A 
vast  amount  of  work  has  recently  been  done  involving  this  method,  largely  by  Haus  and 
his  co-workers.  This  includes  work  on  conditions  necessary  for  self-starting1 14  and  the 
effect  of  third-order  dispersion  on  pulses  in  a  ring  laser.1 15  This  method  has  been  used  to 
develop  a  type  of  laser  known  as  the  stretched  pulse  ring  laser.14  This  method  addressees 
the  observation  that  the  loop  length  must  be  kept  less  than  eight  soliton  periods.  This  was 
done  by  inserting  positive  dispersion  fiber  into  the  loop.  The  pulse  in  this  fiber  will 
spread  out  and  its  phase  shift  per  unit  length  will  fall.  When  the  pulse  reenters  the 
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negative  dispersion  fiber  it  will  be  compressed  again.  These  two  effects  were  balanced 
out  so  that  the  pulse  is  undistorted  in  one  round  trip.  Recently  the  additive  pulse  mode¬ 
locking  formalism  has  been  extended  into  regimes  where  the  pulse  undergoes  larger 
distortions.116  Similar  work  has  been  done  to  analyze  mode-locking  of  a  Ti:sapphire 
laser.117  Other  research  groups  have  examined  conditions  under  which  stable  pulses 
could  be  formed  in  ring  lasers.118,119 

The  side-band  theory  has  been  used  to  study  figure  eight  lasers.  Dennis  and 
Duling  have  made  extensive  investigations  into  this  experimentally.120  They  found  that 
by  reducing  the  product  dispersion  and  the  cavity  length,  they  could  reduce  the  width  of 
the  shortest  pulse  that  could  be  formed  in  their  laser  in  the  same  fashion  as  in  the  ring 
laser.  They  found  that  the  minimum  pulse  formed,  xmjn,  could  be  found  by 

T min  *  ’ 

where  L  is  the  length  of  the  total  cavity  and  p2  is  the  average  dispersion  of  the  cavity. 
The  form  of  this  relationship  can  be  deduced  if  one  inserts  the  expression  of  the 
dispersion  length  into  equation  6.6.  This  yields 


This  theory  concludes  that  the  shortest  pulse  producable  by  a  laser  is  limited  by  the  laser 
length  and  by  its  average  dispersion. 

Sideband  generation  is  of  concern  in  telecommunications  as  well.  The  effect 
occurs  when  solitons  are  maintained  by  a  large  number  of  optical  amplifiers.  It  has  been 
shown  that  sideband  generation  is  reduced  by  introducing  NOLMs  along  with  the 
amplifiers  to  eliminate  the  dispersive  wave  radiation.121  This  implies  that  dispersive 
wave  radiation  should  have  only  a  reduced  role  in  the  lasers  that  will  be  studied  here. 
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It  has  also  been  shown  that  third-order  dispersion  limits  the  length  of  pulses  in 
figure  eight  lasers.  Pulses  could  be  shortened  considerably  by  compensating  for  this 
effect. 


Simulation  Results:  Figure  Eight  Laser 

I  simulated  the  behavior  of  an  NOLM  mode-locked  figure  eight  laser  extensively. 
The  goal  was  to  analyze  the  steady-state  behavior  of  the  laser.  This  was  done  by  inserting 
a  seed  pulse  into  the  laser  shown  in  Figure  6.1.  The  seed  chosen  was  a  hyperbolic  secant 
shaped  pulse  which  provided  maximum  transmutation  in  the  results  given  in  the  last 
chapter.  It  was  assumed  that  this  would  be  a  good  approximation  of  the  steady-state 
pulse  shape.  The  function  of  the  NOLM  was  simulated  in  the  fashion  described  in 
Chapter  5.  The  pulse  was  then  amplified  as  in  Chapter  4.  The  amplifier  was  assumed  to 
have  uniform  gain  along  its  length  and  that  it  fully  recharged  between  passes  of  the  pulse. 
The  effect  of  the  output  coupler  was  simulated  as  a  linear  loss.  For  instance,  for  a  10% 
output  coupler,  the  intensity  of  the  light  at  all  points  was  multiplied  by  0.9.  The  resulting 
pulse  was  then  reintroduced  into  the  NOLM  and  the  cycle  was  continued.  The  pulses 
were  examined  at  the  end  of  the  cycle,  that  is  after  the  output  coupler.  It  was  found  that 
after  fewer  than  1 0  cycles  the  pulses  had  assumed  their  final  shape  in  time  and  frequency 
space.  This  was  determined  by  propagating  a  sample  case  for  10,  20  and  30  cycles  and 
comparing  the  results.  In  general,  the  pulses  were  assumed  to  have  reached  equilibrium  if 
the  total  integrated  energy  in  the  simulation  was  stable  to  within  less  than  5%  for 
successive  cycles. 

To  start  the  investigation  of  the  NOLM  laser,  the  nonbirefringent  case  was 
examined.  The  loop  length  was  assumed  to  be  four  soliton  periods  and  the  directional 
coupler  had  a  60/40  splitting  ratio.  A  10%  output  coupler  was  used  and  the  amplifier  was 
assumed  to  recharge  completely  between  passes  of  the  pulse.  The  effect  of  the 
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polarization  controller  was  ignored.  The  first  fact  that  appears  is  that  one  cannot  use 
simply  any  length  for  the  amplifier.  The  results  in  Figure  6.2,  show  the  maximum  and 
minimum  gain  for  which  stable  pulses  can  be  produced.  In  this  case  gain  is  defined  as 
the  ratio  of  the  power  out  to  the  power  in.  It  is  given  by  e"yL  where  L  is  the  length  of  the 
amplifier  and  y  is  the  gain  parameter.  To  demonstrate  this  pulses  corresponding  to  those 
for  the  maximum  transmission  through  the  NOLM  were  introduced  into  the  laser  and 
allowed  to  pass  through  its  loop  ten  times.  If  the  gain  was  insufficient  the  pulses  would 
decay  away  exponentially.  The  explanation  for  this  behavior  is  simply  that  the  loss  per 
pass  through  the  loop  is  greater  than  the  gain.  When  the  gain  grew  to  be  too  great  a  more 
unexpected  result  occurred.  The  pulse  energy  would  first  oscillate  and  only  decay  away 
at  a  later  time.  An  explanation  for  this  behavior  can  be  derived  from  examining  an 
NOLM  transmission  curve  such  as  Figure  5.2.  If  the  amplification  of  the  pulse  is  too 
great,  the  pulse  energy  will  be  such  that  the  pulse  will  experience  a  great  deal  of  loss  in 
the  NOLM  in  subsequent  passes  and  ultimately  be  extinguished.  Under  these  conditions 
the  pulse  energy  after  each  pass  through  the  system  would  oscillate  several  times  but 
finally  decay  away.  The  amount  of  tolerable  excessive  gain  parameter  appears  to  be 
about  0.1  in  this  case.  This  corresponds  to  an  excess  gain  of  about  10%.  It  should  be 
noted  that  behavior  like  this  is  observed  in  lasers  experimentally.  As  noted  before,  it  has 
been  found  that  figure  eight  lasers  transition  through  several  stability  regimes.  As 
pumping  is  increased,  they  transition  from  having  one  pulse  in  the  loop  to  having  several. 
The  observation  here  of  a  maximum  stable  gain  make  sense  in  this  context.  When  the 
amplification  becomes  too  great,  the  existence  of  this  maximum  will  force  the  laser  to 
produce  two  small  pulses  rather  than  one  pulse  with  greater  energy.  Saturation  due  to  the 
second  pulse  will  eliminate  the  excess  amplification.  When  the  amplification  is  even 
greater,  the  pulse  can  breakup  into  many  pulses,  as  shown  in  the  chapter  on  amplifiers. 
This  would  produce  the  chaotic  bursts  of  pulses  which  are  observed. 

For  short  amplifiers,  the  gain  parameter  needed  corresponds  to  an  amplification  of 
about  10%.  This  means  that  for  these  short  amplifiers,  only  enough  gain  is  needed  to 
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compensate  for  losses  from  the  output  coupler.  For  longer  amplifiers,  there  is  a  loss 
which  is  dependent  on  the  amplifier  length.  Behavior  of  this  sort  has  been  noted  by  a 
number  of  groups  in  studying  lasers  in  a  ring  geometry.  In  that  case,  it  is  thought  to  be 
caused  by  the  repetitive  amplification  of  dispersive  wave  radiation.  The  exact  role  of 
this  phenomenon  in  the  NOLM  is  problematic,  as  noted  above.  Dispersive  wave 
radiation  would  appear  as  a  weak  pedestal  about  the  pulse.  In  the  case  of  the  NOLM 
laser,  however,  this  should  be  eliminated  by  the  switch.  Gordon  also  found  that  loss 
increased  very  rapidly  as  the  loop  length  approached  8  soliton  periods.  My  investigations 
show  that  stable  operation  was  impossible  if  the  amplifier  was  longer  than  about  0.83 
dispersion  lengths. 

The  shape  of  the  pulse  in  the  laser  is  of  interest  in  itself  and  it  also  sheds  light  on 
the  phenomena  described  in  the  last  paragraph.  Figures  6.3  and  6.4  show  how  the  shape 
of  the  pulse  is  altered  as  the  length  of  the  amplifier  is  increased.  In  this  simulation,  the 
pulse  first  passes  through  the  NOLM,  is  amplified,  and  then  is  reduced  by  10%  due  to  the 
output  coupler.  These  pulses  were  the  result  of  10  passes  through  the  laser.  The  pulse  in 
the  Figures  6.3  and  6.4  is  in  the  state  it  has  just  before  reentering  the  NOLM.  The  pulses 
in  both  figures  show  some  pedestal  corresponding  to  some  dispersive  wave  radiation. 

The  pulse  produced  when  the  amplifier  is  0.1  dispersion  length  is  much  closer  to  being  in 
the  the  form 

P0  sech2(x) . 

which  will  be  referred  to  as  a  hyperbolic  secant  shaped  pulse.  The  pulse  for  a  0.8 
dispersion  length  amplifier  much  more  closely  resembles  a  soliton.  In  this  form  the 
argument  of  the  hyperbolic  secant  has  the  form  (t/To)  where  Po  and  To  have  the  standard 
fixed  relationship.  As  we  saw  in  the  last  section,  the  output  pulse  is  almost  the  right 
shape  for  optimal  transmission  by  the  NOLM,  if  it  were  immediately  fed  back  into  it. 

The  transmission  is  very  sensitive  to  pulse  shape.  Hence,  any  pulse  shaping  by  the 
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amplifier  fiber  will  tend  to  create  loss.  A  longer  amplifier  will  obviously  create  more  loss 
and  this  will  create  the  length  dependent  loss  we  observe. 


amplifier  length 


Figure  6.2  Stable  pulse  formation  could  take  place  between  the  curves  in  a 
nonbirefringent  NOLM  laser  with  a  60/40  beamsplitter,  and  a  4  soliton  period  long 
loop  mirror.  The  output  coupling  is  10%.  Amplifier  length  is  given  in  dispersion 

lengths. 


Figure  6.3  Pulse  shape  in  a  nonbirefringent  NOLM  laser  with  a  60/40  beamsplitter, 
and  a  4  soliton  period  long  loop  mirror.  The  output  coupling  is  10%.  The  amplifier 
is  0.1  dispersion  lengths  long.  The  power  of  the  pulse  is  shown  as  a  function  of  time. 

Both  are  in  soliton  units. 
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Figure  6.4  Pulse  shape  in  a  nonbirefringent  NOLM  laser  with  a  60/40  beamsplitter, 
and  a  4  soliton  period  long  loop  mirror.  The  output  coupling  is  10%.  The  amplifier 
is  0.1  dispersion  lengths  long.  The  power  of  the  pulse  is  shown  as  a  function  of 
frequency.  Both  are  in  soliton  units. 
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Figure  6.5  Pulse  shape  in  a  nonbirefringent  NOLM  laser  with  a  60/40  beamsplitter, 
and  a  4  soliton  period  long  loop  mirror.  The  output  coupling  is  10%.  The  amplifier 
is  0.8  dispersion  lengths  long.  The  power  of  the  pulse  is  shown  as  a  function  of  time. 

Both  are  in  soliton  units. 
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Figure  6.6  Pulse  shape  in  a  nonbirefringent  NOLM  laser  with  a  60/40  beamsplitter, 
and  a  4  soliton  period  long  loop  mirror.  The  output  coupling  is  10%.  The  amplifier 
is  0.8  dispersion  lengths  long.  The  power  of  the  pulse  is  shown  as  a  function  of 
frequency.  Both  are  in  soliton  units. 


Simulation  Results:  Dispersion  Balanced  Figure  Eight  Laser 

The  analysis  given  above  suggests  that  to  improve  laser  operation,  one  should 
design  a  laser  in  which  this  pulse  distortion  is  minimized.  One  possible  means  of  doing 


114 


this  is  by  using  a  different  value  of  fh  in  the  amplifier  and  in  the  loop  mirror.  The  value 
of  dispersion  should  be  chosen  so  that  the  output  pulse  from  the  NOLM  is  a  fundamental 
soliton  in  the  amplifier  fiber.  According  to  Equation  2.6,  the  parameters  controlling  the 
relation  between  the  width  of  a  fundamental  soliton  and  its  height  are 
,  YPoT02 


The  parameter  y  is  not  in  general,  easily  changed.  It  is  dependent  on  the  nonlinear 
properties  of  the  fiber  and  on  its  effective  core  area.  The  first  is  a  basic  material  property 
of  silica  and  the  second  tends  to  be  driven  by  the  necessities  of  good  amplifier  design. 
Dispersion,  however,  can  be  controlled  and  dispersion  shifted  fiber  is  available.  It  is 
therefore  plausible  to  consider  changing  this  parameter.  In  the  case  of  a  four  soliton 
period  long  loop,  the  pulse  width  is  not  altered  by  transmission  through  the  NOLM.  This 
means  that  we  want  To  to  be  unaltered  by  propagation  in  the  amplifier  fiber.  The  output 
pulse  has  in  the  case  we  are  examining  a  peak  intensity  2.2  times  that  of  a  unit 
fundamental  soliton.  Since  y  and  To  are  held  constant,  the  only  way  for  Equation  6.7  to 
be  true  in  the  amplifier  is  if  the  dispersion  is  multiplied  by  a  factor  of  about  2.2. 

In  order  to  examine  this  possible  solution,  Equation  2.36,  the  normalized 
nonlinear  Schrodinger  equation  must  be  changed.  The  starting  point  for  this  is  Equation 
2.35,  the  unnormalized  NLSE.  This  has  the  form,  if  we  ignore  higher  order  terms 
dA  a  i  d2A  1  53A  r,  l2  i 

8T+IA+I|,*aF-«f,>aF-  ,y[|A|  AI 

If  the  dispersion  term  is  multiplied  by  a  factor  K  and  the  equation  is  normalized  as  in 
section  2,  the  resulting  equation  is 

5U  K92U  a3U  ,r.  ,2  1 

-l5a^=-N  [|u|  ul'  (6'9) 

In  the  case  we  are  examining,  this  means  that  the  linear  propagation  parameters  will  be 
multiplied  by  the  factor  K. 
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If  propagation  in  the  fiber  amplifier  is  modified  in  this  way,  the  result  is  quite 
striking.  The  stable  gain  region  is  shown  in  Figure  6.5.  This  simulation  did  not  exhibit  a 
maximum  amplifier  length  as  in  Figure  6.1.  This  occurred  even  when  the  simulation  was 
taken  out  to  about  30  times  the  amplifier  length  examined  in  the  noncompensated  case. 

In  addition  this  curve  does  not  exhibit  length  dependent  loss  to  any  great  extent. 


amplifier  length 

Figure  6.7  Stable  pulse  formation  could  take  place  between  the  curves  in  a 
nonbirefringent  dispersion  balanced  NOLM  laser  with  a  60/40  beamsplitter,  and  a  4 
soliton  period  long  loop  mirror.  The  output  coupling  is  10%.  The  dispersion  in  the 
amplifier  is  2.2  times  that  in  the  loop.  Amplifier  length  is  given  in  dispersion 

lengths. 
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The  shape  of  the  pulse  is  shown  in  Figure  6.8  and  6.9.  The  pulse  shape  in  time 
space  has  the  correct  soliton  shape.  This  is  to  be  expected  since  the  amplification  is 
extremely  gradual  and  the  pulse  can  easily  be  shaped  into  a  soliton  by  the  fiber.  The 
pulse  displays  very  little  chirp.  The  pulse  in  frequency  space  is  totally  dominated  by  the 
sidebands.  The  location  of  the  sidebands  is  summarized  in  Figure  6.10.  The  dependent 
variable  is  VF  which  is  defined  to  by 


50 

271 ) 


X 2 


From  Equation  6.5,  we  see  that  the  sideband  theory  predicts  that 


VF  = 


1 


4tc 


—  1 


(6.10) 


(6.11) 


where  m  is  the  order  of  the  side  band.  The  variable  T*  increases  linearly  with  the  order, 
and  Figure  6.10  shows  this  relationship.  The  slopes  of  the  lines  can  be  used  to  compute 
an  effective  length  of  the  cavity,  Leff,  as  the  amplifier  length  is  varied.  This  computation 
is  made  easier  if  the  sidebands  are  compared  to  the  location  of  the  first  sideband.  If  m  is 
the  order  of  the  first  sideband  and  n  is  the  order  of  the  one  to  which  it  is  being  compared, 
one  obtains 


Tf  -  Tf  = 


1  4nmL 


4n2 


-(n-m). 


■'eff 


(6.12) 


By  examining  how  T*  changes  as  the  difference  in  order  changes,  we  may  obtain  an 
effective  cavity  length  in  dispersion  lengths.  We  see  in  figure  6.10  how  the  slope  of  the 
line  changes  as  the  length  of  the  amplifier  changes.  The  effective  length  is  the  length  of 
the  amplifier  that  would  produce  the  sideband  spacing  observered. 

It  is  interesting  to  do  linear  regression  analysis  on  the  relationship  between  the 
effective  length  and  the  length  of  the  amplifier.  This  will  help  us  to  relate  the  effective 
length  of  the  loop  to  the  actual  length.  The  results  of  this  analysis  are  displayed  in  Figure 
6.11.  If  a  linear  regression  is  performed  on  the  data,  one  obtains 
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L«  =2.70L„  -2.62, 


(6.13) 


where  Lamp  is  the  length  of  the  amplifier  in  dispersion  lengths  for  a  soliton  of  unit 
amplitude.  The  relationship  between  the  dispersion  length  and  physical  length  is  the  key 
to  understanding  this  relationship.  From  equation  2.36  we  see  that 

F.JfeL-L. 

0  T2  T 
Ao 

This  indicates  that  we  can  relate  a  change  in  the  peak  intensity  of  a  soliton  directly  to  a 
change  in  its  dispersion  length.  The  peak  intensity  of  the  pulse  is  2.66  for  all  lengths  of 
the  amplifier.  If  the  peak  intensity  after  the  output  coupler  was  2.66,  the  intensity  before 
would  have  been  2.96.  Before  the  output  coupler  this  intensity  would  have  been  2.96. 
This  means  that  the  dispersion  length  of  a  soliton  of  unit  amplitude  should  have  a 
dispersion  length  2.96  times  that  of  the  pulse  at  the  end  of  a  circuit  of  the  laser,  so  the 
pulse  should  have  been  acting  as  if  the  amplifier  was  about  2.96  times  longer  than 
indicated.  This  explains  the  slope  in  equation  6.13.  The  linear  regression  disagrees  with 
this  by  17%.  The  NOLM  is  2tc  dispersion  lengths  long.  This  length  agrees  approximately 
with  the  effective  length  for  short  amplifiers.  This  is  to  be  expected  since  for  a  short 
amplifier  most  of  the  length  of  the  laser  is  in  the  NOLM,  where  the  pulse  has  a  peak 
intensity  of  about  1  in  the  units  used  here. 
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Figure  6.8  Pulse  shape  in  a  nonbirefringent  dispersion  balanced  NOLM  laser  with  a 
60/40  beamsplitter,  and  a  4  soliton  period  long  loop  mirror.  The  output  coupling  is 
10%.  The  dispersion  in  the  amplifier  is  2.2  times  that  in  the  loop.  The  pulse 
intensity  is  given  as  a  function  of  time.  Both  are  in  soliton  units. 
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Figure  6.9  Pulse  shape  in  a  nonbirefringent  dispersion  balanced  NOLM  laser  with  a 
60/40  beamsplitter,  and  a  4  soliton  period  long  loop  mirror.  The  output  coupling  is 
10%.  The  dispersion  in  the  amplifier  is  2.2  times  that  in  the  loop.  Pulse  intensity  is 
given  as  a  function  of  frequency.  Both  are  in  soliton  units. 
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Figure  6.10  Sideband  spacing  in  a  nonbirefringent  dispersion  balanced  NOLM 
laser  with  a  60/40  beamsplitter,  and  a  4  soliton  period  long  loop  mirror.  The  output 
coupling  is  10%.  The  dispersion  in  the  amplifier  is  2.2  times  that  in  the  loop.  The 
parameter  is  defined  in  Equation  6.10.  The  figure  shows  the  effect  of  changing 

the  length  of  the  amplifier. 
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Figure  6.11  Leff  as  computed  by  means  of  Equation  6.12.  The  laser  is  a 
nonbirefringent  dispersion  balanced  NOLM  laser  with  a  60/40  beamsplitter,  and  a  4 
soliton  period  long  loop  mirror.  The  output  coupling  is  10%.  The  dispersion  in  the 

amplifier  is  2.2  times  that  in  the  loop 


The  effect  of  changing  the  splitting  ratio  of  the  directional  coupler  in  the  NOLM 
was  examined.  The  splitting  ratio  in  the  simulation  was  changed  to  55/45.  The  loop  was 
lengthened  to  8  soliton  periods.  This  meant  that  the  input  pulse  energy  for  maximum 
transmission  was  unchanged.  The  two  counter  propagating  pulses  were  approximately 
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unit  solitons.  The  linear  dispersion  was  2.2  times  greater  in  the  amplifier  than  in  the  loop. 
The  resulting  stability  curve  is  shown  in  Figure  6.12.  The  dispersion  balancing  scheme 
works  for  this  case  as  well.  The  laser  will  form  stable  pulses  for  very  long  amplifiers. 


Figure  6.12  Stable  pulse  formation  could  take  place  between  the  curves  in  a 
nonbirefringent  dispersion  balanced  NOLM  laser  with  a  55/45  beamsplitter,  and  a  8 
soliton  period  long  loop  mirror.  The  output  coupling  is  10%.  The  dispersion  in  the 
amplifier  is  2.2  times  that  in  the  loop.  Amplifier  length  is  given  in  dispersion 

lengths. 

The  shape  of  the  pulse  for  this  laser  is  shown  in  Figure  6.13.  This  is  for  a  20 
dispersion  length  amplifier  with  near  minimum  stable  amplification.  There  is  again  very 
little  chirp.  It  should  be  noted  that  side  bands  form  in  the  case  of  the  55/45  directional 
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coupler  in  much  the  same  way  as  they  did  for  the  60/40  directional  coupler.  It  appears 
that  the  sidebands  are  larger  than  those  in  the  case  of  the  laser  with  the  60/40 
beamsplitter.  This  is  of  interest  since  the  counterpropagating  pulses  in  this  case  should 
have  an  initial  form  closer  to  that  of  a  fundamental  soliton.  It  would  seem  that  pulse 
reshaping  in  the  NOLM  is  not  the  major  source  of  dispersive  wave  shedding  of  light  from 
the  pulses  since  the  pulses  were  closer  to  being  solitons  in  this  case. 

The  location  of  the  sidebands  can  also  be  analyzed  as  in  the  case  of  the  60/40 
laser.  The  relationship  between 'T  and  the  sideband  order  is  shown  in  Figure  6.15  for 
lasers  with  various  amplifier  lengths.  The  relationship  between  the  computed  Leff  and  the 
length  of  the  amplifier  is  shown  in  figure  6.16.  The  resulting  formula  for  the  effective 
length,  in  dispersion  lengths,  for  this  laser  is 

Leff=2.29Lamp+ 2.373  (6.14) 

The  peak  intensity  for  this  laser  is  2.52.  Before  the  10%  output  coupler  the  peak  would 
have  been  2.80.  As  was  argued  before  when  describing  equation  6.13,  the  slope  should 
be  around  this  value.  In  this  case  the  difference  between  this  figure  and  the  slope  is  17%. 
The  NOLM  is  now  4n  dispersion  lengths  long  and  this  is  about  the  effective  length  for  a 
laser  with  a  short  amplifier. 
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Figure  6.13  Pulse  shape  in  a  nonbirefringent  dispersion  balanced  NOLM  laser  with 
a  55/45  beamsplitter,  and  a  8  soliton  period  long  loop  mirror.  The  output  coupling 
is  10%.  The  dispersion  in  the  amplifier  is  2.2  times  that  in  the  loop.  Pulse  intensity 
is  given  as  a  function  of  time.  Both  are  in  soliton  units. 
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Figure  6.14  Pulse  shape  in  a  nonbirefringent  dispersion  balanced  NOLM  laser  with 
a  55/45  beamsplitter,  and  a  8  soliton  period  long  loop  mirror.  The  output  coupling 
is  10%.  The  dispersion  in  the  amplifier  is  2.2  times  that  in  the  loop.  Pulse  intensity 
is  given  as  a  function  of  frequency.  Both  are  in  soliton  units. 
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Figure  6.15  Sideband  spacing  in  a  nonbirefringent  dispersion  balanced  NOLM 
laser  with  a  55/45  beamsplitter,  and  a  8  soliton  period  long  loop  mirror.  The  output 
coupling  is  10%.  The  dispersion  in  the  amplifier  is  2.2  times  that  in  the  loop.  The 
parameter  'F  is  defined  in  Equation  6.10.  The  figure  shows  the  effect  of  changing 

the  amplifier  length. 
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Figure  6.16  Leff  as  computed  by  means  of  Equation  6.12.  The  laser  is  a 
nonbirefringent  dispersion  balanced  NOLM  laser  with  a  55/45  beamsplitter,  and  a  8 
soliton  period  long  loop  mirror.  The  output  coupling  is  10%.  The  dispersion  in  the 
amplifier  is  2.2  times  that  in  the  loop.  All  lengths  are  in  dispersion  lengths. 
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Figure  6.17  Stable  pulse  formation  could  take  place  between  the  curves  in  a 
nonbirefringent  dispersion  balanced  NOLM  laser  with  a  60/40  beamsplitter,  and  a  4 
soliton  period  long  loop  mirror.  The  output  coupling  is  30%.  The  dispersion  in  the 
amplifier  is  3.0  times  that  in  the  loop.  Amplifier  length  is  given  in  dispersion 

lengths. 


The  effect  of  increasing  the  output  coupling  was  examined  by  changing  the 
coupling  to  30%.  The  central  beamsplitter  ratio  was  left  at  60/40.  The  loop  was  4  soliton 
periods  long.  Finding  a  value  of  birefringence  for  the  amplifier  that  would  work  in  this 
case  proved  to  be  challenging.  The  laser  was  first  simulated  with  a  dispersion 
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compensation  value  of  2.2.  This  was  increased  and  the  rate  at  which  the  pulse  in  the  laser 
was  amplified  was  observed.  The  amplification  was  highest  when  the  dispersion  in  the 
amplifier  was  three  times  that  in  the  loop.  The  result  of  simulating  a  laser  with  this 
configuration  is  shown  in  Figure  6.17.  Again,  long  cavities  can  be  achieved.  The  need  to 
increase  the  dispersion  in  this  configuration  is  to  be  expected.  The  amplifier  section  of 
the  laser  is  immediately  followed  by  the  output  coupler.  Since  the  pulse  must  be 
amplified  enough  so  that  it  has  sufficient  energy  to  be  transmitted  by  the  NOLM  even 
after  losing  30%  of  its  energy  to  the  output  coupler.  This  means  that  a  value  of  dispersion 
must  be  chosen  such  as  to  obtain  the  right  pulse  shape  even  after  amplification. 

The  sidebands  in  this  case  are  extremely  pronounced  as  can  be  seen  in  Figure 
6.19.  This  occurs  due  to  the  larger  perturbations  the  laser  imposes  on  the  soliton  like 
pulses.  This  make  the  sideband  analysis  easy  and  the  sideband  spacing  information  is 
summarized  in  Figure  6.20.  Linear  regression  was  used  to  find  the  effective  length  as  a 
function  of  the  amplifier  length.  This  process  yields  the  curve  in  Figure  6.21 .  Again,  all 
these  functions  appear  to  be  approximately  linear.  The  form  of  the  line  in  Figure  6.21  is 

Leff=3.29Lamp- 1.1087.  (6.15) 

The  pulse  peak  intensity  is  2.60.  This  is  just  after  the  output  coupler  and  before  the  30% 
output,  so  before  the  couple  the  intensity  would  have  been  3.71 .  This  means  that  the 
nonlinear  length  should  have  been  the  same.  This  value  is  about  19%  larger  than  the 
slope  in  Equation  6.14.  It  is  interesting  to  note  that  the  slopes  were  all  smaller  than  the 
expected  values  based  on  the  peak  laser  intensities  by  about  the  same  amount.  This  may 
imply  that  there  is  in  fact  some  simple  explanation  possible  for  the  disagreement. 

As  was  mentioned  before,  the  length  of  fiber  lasers  has  been  assumed  to  be 
limited  to  about  8  soliton  periods.  Using  this  length  in  Equation  6.15,  we  see  that  this 
corresponds  to  an  effective  length  of  4.02  dispersion  lengths.  When  Figure  6.17  is 
examined,  this  length  is  found  to  be  where  a  large  broadening  of  the  stable  gain  region 
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occurs.  Equations  6.13  and  6.15  can  be  used  to  compute  the  amplifier  lengths 
corresponding  to  8  soliton  periods  in  the  cases  seen  in  Figures  6.7  and  6.12.  The 
computed  values  are  5.3  and  5.2  dispersion  lengths  respectively.  In  both  cases,  the  lasers 
appear  to  have  some  change  in  the  stable  gain  curve  near  this  point  though  its  character  is 
much  less  dramatic  than  in  the  laser  with  high  output.  This  indicates  that  while  these 
lasers  do  not  cease  functioning  with  an  effective  length  of  8  soliton  periods  they  do 
undergo  some  shift  in  their  operation. 


time 

- pulse  shape  —  — phase 

Figure  6.18  Pulse  shape  in  a  nonbirefringent  dispersion  balanced  NOLM  laser  with 
a  60/40  beamsplitter,  and  a  4  soliton  period  long  loop  mirror.  The  output  coupling 
is  30%.  The  dispersion  in  the  amplifier  is  3.0  times  that  in  the  loop.  Pulse  shape  is 
given  as  a  function  of  time.  Both  are  in  soliton  units. 
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Figure  6.19  Pulse  shape  in  a  nonbirefringent  dispersion  balanced  NOLM  laser  with 
a  60/40  beamsplitter,  and  a  4  soliton  period  long  loop  mirror.  The  output  coupling 
is  30%.  The  dispersion  in  the  amplifier  is  3.0  times  that  in  the  loop.  Pulse  shape  is 
given  as  a  function  of  frequency.  Both  are  given  in  soliton  units. 
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Figure  6.20  Sideband  spacing  in  a  nonbirefringent  dispersion  balanced  NOLM 
laser  with  a  60/40  beamsplitter,  and  a  4  soliton  period  long  loop  mirror.  The  output 
coupling  is  30%.  The  dispersion  in  the  amplifier  is  3.0  times  that  in  the  loop.  The 
parameter  'F  is  defined  in  Equation  6.10.  The  figure  shows  the  effect  of  changing 

the  amplifier  length. 
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amplifier  length 


Figure  6.21  Leff  as  computed  by  means  of  Equation  6.12.  The  laser  is  a 
nonbirefringent  dispersion  balanced  NOLM  laser  with  a  60/40  beamsplitter,  and  a  4 
soliton  period  long  loop  mirror.  The  output  coupling  is  30%.  The  dispersion  in  the 
amplifier  is  3.0  times  that  in  the  loop.  Lengths  are  given  in  dispersion  lengths. 
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Two  situations  were  simulated  in  which  the  lasers  were  assumed  to  have 
birefringent  fiber.  The  first  was  the  case  in  which  all  of  the  starting  pulse  energy  was  on 
the  slow  birefringent  axis  and  the  plane  of  the  loop  tilted  by  tt/4  radians  from  this 
direction.  In  the  second,  the  pulse  was  a  vector  soliton  and  had  equal  amplitudes  along 
both  axes.  The  plane  of  the  loop  was  along  the  slow  axis.  These  cases  were  chosen 
because  they  appeared  to  offer  the  possibility  of  controlling  the  NOLM  switch  by 
adjusting  the  polarization  controller.  This  was  discussed  in  the  chapter  on  the  NOLM. 

The  results  are  displayed  in  Figures  6.22-6.24.  The  dispersion  in  the  amplifier 
was  made  to  be  2.2  times  that  in  the  loop.  The  pulse  in  the  laser  is  again  a  nearly 
unchirped  soliton.  The  shift  of  the  pulse  center  in  time  is  due  to  the  fact  that  it  is 
propagating  predominately  on  the  slow  axis.  Essentially  none  of  the  light  leaked  into  the 
orthogonal  polarization.  In  frequency  space,  the  location  of  the  maximum  of  the  pulse 
has  not  changed,  but  the  pulse  has  become  asymmetric  and  the  location  of  the  sidebands 
have  shifted. 

A  laser  with  a  vector  soliton  trial  pulse  was  simulated.  The  results  are  in  Figures 
6.25-6.27.  The  dispersion  in  the  amplifier  was  again  increased  by  2.2.  The  laser  pulse  is 
found  to  maintain  its  polarization.  In  frequency  space,  one  observes  that  the  two 
polarizations  couple  with  each  other.  This  is  also  the  explanation  of  the  chirp  that 
develops  in  the  phase  and  the  shift  observed  in  the  central  frequencies. 

When  a  phase  shift  of  7t/4  was  imposed  with  the  polarization  controller  no  stable 
operating  configuration  could  be  found.  This  was  true  for  both  of  the  cases  examined  in 
which  birefringence  was  considered.  As  described  in  the  last  section,  the  polarization 
controller  reshapes  the  pulse.  The  same  effects  that  lead  to  poor  two-pass  transmission  in 
the  case  of  the  NOLM  lead  to  excessively  high  loss  in  the  laser.  It  was  found  that  this 
loss  could  not  be  compensated  with  additional  amplification.  The  pulse  would  stabilize 
for  a  pass  and  then  decay  away. 
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Figure  6.22  Stable  pulse  formation  could  take  place  between  the  curves  in  a 
birefringent  NOLM  laser  with  a  60/40  beamsplitter,  and  a  4  soliton  period  long  loop 
mirror.  The  output  coupling  is  10%.  The  birefringence  parameter  is  8=0.1.  The 
input  pulse  and  the  slow  axis  of  the  polarization  controller  are  aligned  with  the  slow 
birefringent  axis.  The  plane  of  the  loop  has  a  -tt/4  angle  with  the  slow  axis.  The 
dispersion  in  the  amplifier  is  2.2  times  that  in  the  loop.  The  amplifier  length  is  given 

in  dispersion  lengths. 
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Figure  6.23  Pulse  shape  in  a  birefringent  NOLM  laser  with  a  60/40  beamsplitter, 
and  a  4  soliton  period  long  loop  mirror.  The  output  coupling  is  10%.  The 
birefringence  parameter  is  8=0.1.  The  input  pulse  and  the  slow  axis  of  the 
polarization  controller  are  aligned  with  the  slow  birefringent  axis.  The  plane  of  the 
loop  has  a  -tc/4  angle  with  the  slow  axis.  The  dispersion  in  the  amplifier  is  2.2  times 
that  in  the  loop.  Pulse  intensity  is  given  as  a  function  of  time.  Both  are  given  in 

soliton  units. 
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Figure  6.24  Pulse  shape  in  a  birefringent  NOLM  laser  with  a  60/40  beamsplitter, 
and  a  4  soliton  period  long  loop  mirror.  The  output  coupling  is  10%.  The 
birefringence  parameter  is  5=0.1.  The  input  pulse  and  the  slow  axis  of  the 
polarization  controller  are  aligned  with  the  slow  birefringent  axis.  The  plane  of  the 
loop  has  a  -7t/4  angle  with  the  slow  axis.  The  dispersion  in  the  amplifier  is  2.2  times 
that  in  the  loop.  Pulse  intensity  is  given  as  a  function  of  frequency.  Both  are  in 

soliton  units. 
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Figure  6.25  Stable  pulse  formation  could  take  place  between  the  curves  in  a 
birefringent  NOLM  laser  with  a  60/40  beamsplitter,  and  a  4  soliton  period  long  loop 
mirror.  The  output  coupling  is  10%.  The  birefringence  parameter  is  8=0.1.  The 
input  pulse  and  the  slow  axis  of  the  polarization  controller  are  at  an  angle  of  tx/4. 
with  the  slow  birefringent  axis.  The  plane  of  the  loop  is  aligned  with  the  slow  axis. 
The  dispersion  in  the  amplifier  is  2.2  times  that  in  the  loop 
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Figure  6.26  Pulse  shape  in  a  birefringent  NOLM  laser  with  a  60/40  beamsplitter, 
and  a  4  soliton  period  long  loop  mirror.  The  output  coupling  is  10%.  The 
birefringence  parameter  is  8=0.1.  The  input  pulse  and  the  slow  axis  of  the 
polarization  controller  are  at  an  angle  of  rc/4.  with  the  slow  birefringent  axis.  The 
plane  of  the  loop  is  aligned  with  the  slow  axis.  The  dispersion  in  the  amplifier  is  2.2 
times  that  in  the  loop.  Pulse  intensity  is  given  as  a  function  of  time.  Both  are  given 

in  soliton  units. 
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Figure  6.27  Pulse  shape  in  a  birefringent  NOLM  laser  with  a  60/40  beamsplitter, 
and  a  4  soliton  period  long  loop  mirror.  The  output  coupling  is  10%.  The 
birefringence  parameter  is  5=0.1.  The  input  pulse  and  the  slow  axis  of  the 
polarization  controller  are  at  an  angle  of  7i/4.  with  the  slow  birefringent  axis.  The 
plane  of  the  loop  is  aligned  with  the  slow  axis.  The  dispersion  in  the  amplifier  is  2.2 
times  that  in  the  loop.  Pulse  intensity  is  given  as  a  function  of  frequency.  Both  are 

given  in  soliton  units. 


The  conclusion  of  the  results  appears  to  be  that  the  NOLM  figure  eight  laser  can 
function  in  a  number  of  ways  depending  upon  the  amount  of  loss  in  the  system.  The  loss 
can  be  introduced  into  the  system  in  a  number  of  ways.  In  the  results  above,  we  have 
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seen  that  it  can  be  introduced  by  pulse  reshaping  in  the  amplifier,  by  birefringence 
reducing  transmission  through  the  NOLM,  or  by  linear  loss  in  the  output  coupler.  If  loss 
in  all  of  these  are  low,  then  it  is  possible  to  make  the  laser  cavity  very  long.  Sidebands 
will  build  up,  but  they  will  not  be  great  enough  to  make  the  laser  stop  working.  If  the 
laser  has  more  loss  in  it,  then  it  will  function  in  a  fashion  described  before  by  other  in 
which  the  growth  of  the  side  bands  creates  a  limit  to  how  long  the  laser  cavity  can  be 
made. 

Haus116  introduced  the  idea  of  inserting  normal  dispersion.  This  was  described  as 
a  means  of  reshaping  the  pulse  and  they  said  that  they  optimized  the  results  by  computer 
simulation.  Dulling  did  extensive  experiments  based  on  the  idea  of  reducing  total 
dispersion  in  the  cavity.  The  new  idea  which  is  introduced  here  is  that  in  a  figure  eight 
laser  this  balancing  can  be  done  by  actually  increasing  dispersion.  Dennis  andDulling120 
did  not  invoke  the  idea  of  pulse  reshaping  nor  did  he  describe  where  he  inserted  the 
positive  dispersion  fiber.  I  have  been  able  to  take  the  pulse  reshaping  ideas  of  Haus  and 
extend  them  in  a  creative  way  to  deal  with  the  specific  dynamics  of  the  figure  eight 
configuration. 

The  most  straightforward  experimental  comparison  is  to  the  work  of  Dennis  and 
Dulling.120  They  did  experiments  with  an  NALM  mode-locked  figure  eight  laser.  They 
experimentally  confirmed  a  relationship  between  the  total  length  of  the  laser  and  the 
minimum  producable  pulse  width.  The  relationship  is  given  by 

^min  00  , 

as  above.  This  result  is  consistent  with  the  observation  here  that  there  is  a  limit  to  how 
long  the  amplifier  can  be  made.  Since  the  length  is  given  in  soliton  units,  the  length  of 
the  cavity,  and  the  width  of  the  pulses  is  linked  by  this  equation.  They  also  observed  the 
formation  of  sidebands.  The  sideband  theory  has  been  verified  for  ring  lasers. 111  It 
appears  that  the  predictions  of  the  simulations  given  here  agree  as  long  as  the  effective 
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length  of  the  loop  is  adjusted  to  take  into  account  the  fact  that  the  pulse  has  a  different 
energy  in  different  parts  of  the  loop. 


143 


Chapter  Seven  :  Conclusions 

The  research  described  in  this  thesis  has  extended  understanding  of  fiber  optic 
devices  and  lasers.  The  work  has  been  principally  concerned  with  applying  existing 
models  to  give  a  more  detailed  understanding  of  their  operation.  Much  of  this  work  has 
centered  upon  birefringence  in  optical  fibers.  Some  important  new  ideas  have  emerged 
about  nonbirefringent  figure  eight  lasers,  as  well. 

New  results  have  been  obtained  concerning  the  amplification  of  pulses  in 
birefringent  optical  amplifiers.  The  simulations  showed  that,  for  high  amplification,  the 
pulse  broke  up  into  many  pulses,  and  that  the  pulses  formed  tended  to  couple  with  pulses 
formed  in  the  orthogonal  polarization.  This  results  in  a  train  of  coupled  solitary  waves 
propagating  together.  This  coupling  can  be  broken  by  a  large  amplifier  saturation 
parameter. 

The  simulation  of  the  nonlinear  optical  loop  mirror  (NOLM)  was  also  extended 
into  the  birefringent  case.  Specifically,  I  studied  the  effect  of  birefringence  on  the 
function  of  the  NOLM.  Many  of  the  same  phenomena  were  found  as  in  the 
nonbirefringent  case.  It  was  found  that  the  length  of  the  output  pulse  relative  to  the  input 
pulse  tends  to  increase  with  the  length  of  the  loop;  however,  the  energy  required  for 
transmission  decreases.  I  examined  the  effect  of  adjusting  the  polarization  controller  in 
the  loop.  The  controller  reduces  the  energy  required  for  switching,  and  increases  low 
intensity  transmission.  Unfortunately,  it  tends  to  distort  pulses  as  well.  This  reduces  the 
transmission  of  the  pulses  when  reintroduced  into  a  controller,  and  hence,  the  utility  of 
the  device.  The  problem  is  reduced  if  the  birefringence  in  the  loop  is  small.  Since  the 
values  of  birefringence  examined  here  are  fairly  high,  the  technique  may  still  have  utility. 

The  work  on  the  amplifier  and  on  the  NOLM  was  combined  to  study  a  figure 
eight  laser  mode-locked  by  an  NOLM.  The  eight  soliton  period  length  limit  found  by 
other  groups  experimentally,  and  analytically  by  other  means,  reemerged  in  this 
simulation.  A  maximum  gain  for  a  stable  pulse  also  was  found.  This  limit  puts  a  bound 


144 


on  how  powerful  a  pulse  can  become  in  a  laser  of  this  type.  It  also  explains  why  pulses 
always  break  up  into  many  smaller  pulses,  when  pump  power  is  increased. 

It  was  demonstrated,  for  the  first  time,  that  the  limit  of  the  length  of  the  loop  could 
be  eliminated  by  balancing  the  dispersion  in  the  two  parts  of  the  figure  eight  laser.  If  the 
dispersion  in  the  amplifier  section  is  increased,  so  that  the  pulse  is  nearly  a  soliton  in  both 
parts  of  the  laser,  then  the  cavity  can  be  made  as  long  as  one  wishes.  This  method 
eliminates  a  design  constraint  for  this  type  of  laser. 

The  pulses  in  this  simulation  displayed  frequency  side-bands  of  the  kind  that  are 
experimentally  found.  Their  location  agreed  with  experimental  observations  and 
analytical  predictions.  An  effective  length  of  the  cavity  was  predicted.  This  effective 
length  could  be  related  to  the  physical  length  of  the  cavity,  and  to  the  relative  pulse 
intensities  in  the  different  parts  of  the  laser. 

A  laser  with  birefringence  was  simulated.  It  could  be  made  to  function  when  the 
polarization  controller  applied  no  phase  delay  along  one  birefringent  axis  relative  to  the 
other.  When  a  phase  delay  was  imposed,  however,  the  controller  created  enough  pulse 
distortion  to  cause  the  laser  to  never  produce  stable  pulses.  This  implies  that  the  behavior 
of  the  NOLM  laser  can  not  be  controlled  in  this  manner,  unless  the  birefringence  of  the 
fiber  is  very  small.  If  the  birefringence  is  large,  the  controller  must  be  adjusted  so  as  to 
eliminate  this  phase  delay. 

The  idea  of  dispersion  balanced  figure  eight  lasers  has  not  been  examined  at  all  in 
the  past.  The  method  of  increasing  dispersion  in  one  part  of  the  fiber  actually  goes  in  the 
opposite  direction  of  the  methods  used  by  other  researchers.  At  the  present  time 
experiments  to  verify  the  predictions  of  the  model  are  being  considered. 
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